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Abstract 



An introductory overview on Light-Front quantization, with some empha- 
sis on recent achievements, is given. Light-Front quantization is the most 
promising and physical tool to study deep inelastic scattering on the basis of 
quark gluon degrees of freedom. The simplified vacuum structure (nontriv- 
ial vacuum effects can only appear in zero-mode degrees of freedom) and the 
physical basis allows for a description of hadrons that stays close to intuition. 
Recent progress has ben made in understanding the connection between ef- 
fective LF Hamiltonians and nontrivial vacuum condesates. Discrete Light- 
Cone Quantization, the transverse lattice and Light-Front Tamm-Dancoff (in 
combination with rcnormalization group techniques) are the main tools for 
exploring LF-Hamiltonians nonperturbatively. 
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Chapter 1 
Introduction 



1.1 Deep Inelastic Scattering 

More than twenty years after quantum chromodynamics (QCD) was intro- 
duced as a microscopic theory of strong interactions, very httle is known 
about its solutions. At least in principle, it should be possible to describe the 
interaction of nucleons with external probes using quark and gluon degrees 
of freedom on the basis of QCD. So far, however, the extreme complexity of 
this theory has slowed any progress in this direction considerably. 

Deep inelastic lepton-nucleon scattering (DIS) provides access to quark 
and gluon degrees of freedom in nucleons and nuclei. In these experiments 
one shoots high energy leptons (e.g. electrons) at a hadronic target (usu- 
ally protons or nuclei) and measures the energy and momentum transfer to 



the target by detecting the final state lepton (Fig. LI). The hadronic final 
state X is not measured (usually the nucleon is destroyed in these reactions 
and the hadronic final state consists of many particles). Because of the ex- 
tremely large momentum transfer to the target (typical momentum transfers 
in DIS experiments are several GeV/c or more), the inclusive cross sections 
are dominated by single particle response functions along the light-cone. To 
illustrate this let us use the optical theorem which relates the differential 
lepton nucleon cross section to the imaginary part of the forward Compton 
amplitude f^n 83|| (Fig. |1.2|) . One finds 



dndE' \E J 2n ^ ' ' 

where E, E' are the energies of the initial and final lepton. q = k — k' 
is the four momentum transfer of the lepton on the target and 
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Figure 1.1: Inclusive process e + — > e ' + X, where X is an unidentified 
hadronic state. 

Ifj^iy = 2k^kl + 2kyk'^ + q^g^u is the leptonic tensor. The hadronic tensor 

q) = ^T.J ^^^'■^(^' ('^.(^)^^(O)) \P, S) (1.2) 

{S is the spin of the target proton) contains all the information about the 
parton substructure of the target proton. 

In the Bjorken limit (Q^ = —q"^ oo, P ■ q —>■ oo, Xbj = Q^/2P ■ q 
fixed), deep inelastic structure functions exhibit Bjorken scaling: up to 
kinematical coefficients, the hadronic tensor ( |1.2|) depends only on xbj 
but no longer on (within perturbatively calculable logarithmic correc- 
tions). In order to understand this result, it is convenient to introduce light- 
front variables a^: = = (a° ± a^) / so that the scalar product reads 
a-b = a^b~^ + a„6^ — aj_bj_ = a^b_ + a_6+ — aj_b±. Furthermore let us choose 
a frame where q± = 0. The Bjorken limit corresponds to and g_ fixed, 
while g+ — >■ oo. Bjorken scaling is equivalent to the statement that the struc- 
ture functions become independent of g+ in this limit (again up to trivial 
kinematic coefficients). In this limit, the integrand in Eq.( |1.2D contains the 
rapidly oscillating factor exp{iq^x^) , which kills all contributions to the in- 



tegral except those where the integrand is singular ||Ja 85|| . Due to causality, 
the integrand must vanish for x"^ = 2x~^x~ — < and the current product 
is singular at x~^ = 0, x± = 0. The leading singularity can be obtained from 
the operator product expansion by contracting two fermion operators in the 
product T {Jfj^{x)J,y{0)) = T (^■?/'(x)7^'?/'(a;)?/'(0)7,^'?/'(0)^, yielding a nonlocal 
term bilinear in the fermion field multiplying a free (asymptotic freedom!) 
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Figure 1.2: Inclusive lepton nucleon cross section expressed in terms of the 
imaginary part of the forward Compton amphtude. For = —q^ — >■ oo 
only the 'handbag diagram' (both photons couple to the same quark) sur- 
vives. The 'crossed diagram' (the two photons couple to different quarks) is 
suppressed because of wavefunction effects. 

fermion propagator from to a; which gives rise to the abovementioned sin- 
gularity structure ||CL 84|| . The = x± = Q dominance in the integral has 
two consequences. First it explains Bjorken scaling, because g+ enters the 
hadronic tensor only via the term x^q+ in the exponent and for x+ = the 
g+ dependence drops out. Second, and this is very important for practical 
calculations, the parton distributions, i.e. the Bjorken scaled structure func- 
tions, can be expressed in terms of correlation functions along the light-front 
space direction x~ . For example, for the spin averaged parton distribution 
one obtains 

/dT~ - 
^(P|V'(0)7-V'(a;")|P) exp(?P_a;-XB,), (1.3) 

The physical origin of this result can be understood as follows. Consider 
again the virtual forward Compton amplitude (Fig. |1.2|) . In principle, the 
photons in the first and second interaction in Fig. |1.2| can couple to the same 
as well as to different quarks in the target. However, the hadronic wave- 
function can only absorb momenta which are of the order of the QCD-scale 
{^QCD ~ 200MeV"). Therefore, in the limit of large momentum transfer, 
only such diagrams survive where the two photons in Fig. [T]^ couple to 
the same quark. All other diagrams have large momenta flowing through 
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the wavefunction or they involve extra hard gluon exchanges which results 
in their suppression at large Q^. The large momentum transfer is also im- 
portant because of asymptotic freedom. Since as'(Q^) ~ 1/ log (^Q^/Ag^^^ , 
the running coupling constant of QCD, goes to zero for large Q^, all inter- 
actions of the struck quark can be neglected and it propagates essentially 
without interaction between the two photon-vertices. Furthermore, since the 
momentum transfer is much larger than the masses of the quarks in the tar- 
get, the struck quarks propagation between becomes ultra-relativistic, i.e. it 
moves exceedingly close to the light cone = 0. Due to the high-energy 
nature of the scattering, the relativistic structure function is a LF correlation 
|Ja 72 , Ji 93 1 . Already at this point it should be clear that LF-coordinates 



play a distinguished role in the analysis of DIS experiments — a point which 
will become much more obvious after we have introduced some of the formal 
ideas of LF quantization. 



1.2 Advantages of Light-Front Coordinates 

LF quantization is very similar to canonical equal time (ET) quantization 
[pi 49|| (here we closely follow Ref. |Po 91|| ). Both are Hamiltonian formu- 



lations of field theory, where one specifies the fields on a particular initial 
surface. The evolution of the fields off the initial surface is determined by 
the Lagrangian equations of motion. The main difference is the choice of 
the initial surface, = for ET and a;"^ = for the LF respectively. In 
both frameworks states are expanded in terms of fields (and their derivatives) 
on this surface. Therefore, the same physical state may have very different 
wavefunction^ in the ET and LF approaches because fields at a;° = provide 
a different basis for expanding a state than fields at x"*" = 0. The reason is 
that the microscopic degrees of freedom — field amplitudes at x° = versus 
field amplitudes at x"^ = — are in general quite different from each other 
in the two formalisms. 

This has important consequences for the practical calculation of parton 
distributions ( |1.3|) which are real time response functions in the equal time 
formalism. [| In order to evaluate Eq.( |1.3D one needs to know not only 
the ground state wavefunction of the target, but also matrix elements to 
excited states. In contrast, in the framework of LF quantization, parton 
distributions are correlation functions at equal LF-time , i.e. within the 
initial surface x"*" = and can thus be expressed directly in terms of ground 
state wavefunctions (As a reminder: ET wavefunctions and LF wavefunctions 

^By "wavefunction" we mean here the collection of all Fock space amplitudes. 
^The arguments of and if) in Eq.(1.3) have different time components! 



6 



are in general different objects). In the LF framework, parton distributions 
f{xBj) can be easily calculated and have a very simple physical interpretation 
as single particle momentum densities, where xbj measures the fraction of 
momentum carried by the hadron f\ 

parton 

^Bj = phadron " iX''^) 

Although DIS is probably the most prominent example for practical ap- 
plications of LF coordinates, they prove useful in many other places as 
well. For example, LF coordinates have been used in the context current 
algebra sum rules in particle physics ||FF 65|| . Another prominent example 



is form factors, where moments of the wave function along the LF deter- 
mine the asymptotic falloff at large momentum transfer ||BL 80|| . More re- 



cently, LF quantization found applications in inclusive decays of heavy quarks 
[ |BD+ 92i PF92Hi ^^Z+ 95|| . 



From the purely theoretical point of view, various advantages of LF quan- 
tization derive from properties of the ten generators of the Poincare group 
(translations P^, rotations L and boosts K) |pi 49| , [Ho 91|| . Those generators 



which leave the initial surface invariant (P and L for ET and P_, Pi, L3 and 
K for LF) are "simple" in the sense that they have very simple representa- 
tions in terms of the fields (typically just sums of single particle operators). 
The other generators, which include the "Hamiltonians" (Pq, which is con- 
jugate to x° in ET and P_|_, which is conjugate to the LF-time x"*" in LF 
quantization) contain interactions among the fields and are typically very 
complicated. Generators which leave the initial surface invariant are also 
called kinematic generators, while the others are called dynamic generators. 
Obviously it is advantageous to have as many of the ten generators kinematic 
as possible. There are seven kinematic generators on the LF but only six in 
ET quantization. 

The fact that P_, the generator of x~ translations, is kinematic (obvi- 
ously it leaves x"^ = invariant!) and positive has striking consequences for 
the LF vacuum ||Ho 91|| . For free fields = tv? implies for the LF energy 
p+ = (m^ + pj_) /2p_. Hence positive energy excitations have positive p_. 
After the usual reinterpretation of the negative energy states this implies 
that p_ for a single particle is positive (which makes sense, considering that 
p_ = (po ~ Ps) / "\/2)- P- being kinematic means that it is given by the sum 

■^In DIS with nonrelativistic kinematics (e.g. thermal neutron scattering off hquid 
'^He) one also observes scaling and the structure functions can be expressed in terms of 
single particle response functions. However, due to the different kinematics, nonrelativistic 
structure functions at large momentum transfer are dominated by Fourier transforms of 
equal time response functions, i.e. ordinary momentum distributions. 



7 



of single particle p_ . Combined with the positivity of p_ this implies that the 
Fock vacuum (no particle excitations) is the unique state with P_ = 0. All 
other states have positive P_. Hence, even in the presence of interactions, 
the LF Fock vacuum does not mix with any other state and is therefore an 
exact eigenstate of the LF Hamiltonian P+ (which commutes with If 
one further assumes parity invariance of the ground state this implies that 
the Fock vacuum must be the exact ground state of the fully interacting 
LF quantum field theory. Q In sharp contrast to other formulations of field 
theory, the LF- vacuum is trivial! This implies a tremendous technical advan- 
tage but also raises the question whether nonperturbative LF-field theory is 
equivalent to conventional field theory, where nonperturbative effects usually 
result in a highly nontrivial vacuum structure. This very deep issue will be 
discussed in more detail in Chapter |[ 

Dirac was the first who had the idea to formulate field theory in LF- 
coordinates ||Di 49|| .p| In this remarkable work (almost 20 years before scaling 
was discovered in deep inelastic lepton nucleon scattering !) he has shown 
that it should in principle be possible to formulate a consistent quantum 
theory on the LF. This work laid the basis for all further developments, but 
left many details open. The main issues are the structure of the vacuum, 
renormalization and practical algorithms for solution. 



1.3 Outline 

There are many similarities between the formal steps in ET quantization and 
LF quantization. In Chapter ^ we will explain the basic steps in construct- 
ing LF Hamiltonians and give examples for scalar fields, fermions and gauge 
fields. The vacuum on the LF is very controversial. On the one hand simple 
kinematical arguments seem to show that in LF field theory the vacuum of 
interacting field theories is the same as the free field theory vacuum (all inter- 
actions turned off). In QCD we know that chiral symmetry is spontaneously 
broken. It is up to now unclear whether a LF Hamiltonian, with its trivial 
vacuum, is capable of describing this physics. We will elaborate on this point 
in Chapter ^ Renormalization is an issue because the LF-approach to field 
theory is not manifestly covariant. Thus UV-divergences (which occur on 
the LF as they do everywhere in quantum field theory) are not necessarily 

^Practical calculations show that typical LF Hamiltonians are either unbounded from 
below or their ground state is indeed the Fock vacuum. 

^Later, a similar framework was developed independently on the basis of a Lorentz 



frame ("the infinite momentum frame") that moves with v c [FF 65, 3u 68 



KS 70, BK+ 71, BR+ 73 



We 6C 
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the same for all Lorentz components of a particular operator under consider- 
ation. Clearly this requires a more complex counterterm structure to render 
the theory finite and to restore Lorentz invariance for physical observables 
(see Chapter^). Despite certain technical simplifications, field theory on the 
LF is a priori still an enormously complex many body problem. In particular 
in QCD one knows from DIS experiments that the nucleon consists not only 
of the three valence quarks, but that sea quark pairs and gluons are a sig- 
nificant, if not dominant, component of the nucleon's LF wavefunction, i.e. 
one should not expect that the LF wavefunctions of ground state hadrons in 
QCD are simple. Recent attempts to cast LF bound state problems into a 
form that can be solved on a computer will be described in Chapter ^. 
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Chapter 2 

Canonical Quantization 



2.1 Quantization in Light-Front Coordinates 

In this chapter, the formal steps for quantization on the hght-front are pre- 
sented. For pedagogical reasons this will be done by comparing with conven- 
tional quantization (with as "time"). On the one hand this shows that 
the basic steps in the quantization procedure in LF and in ET formalism 
are in fact very similar. More importantly, however, we will thus be able to 
highlight the essential differences between these two approaches to quantum 
field theory more easily. 

In the context of canonical quantization one usually starts from the action 

S = J d^xC. (2.1) 

[C = C{(f), d^(j))) After selecting a time direction r Q one forms the momenta 
which are canonically conjugate to (f) 

and postulates canonical commutation relations between fields and corre- 
sponding momenta at equal "time" r (Table [2.11 ). f\ 

In the next step one constructs the Hamilton operator and the other 
components of the momentum vector. Thus one has completely specified the 

^Here r may stand for ordinary time as well as for LF time x+ = (a;° + .t"^) /V2 or 
any other (not space-like) direction. 

^The canonical quantization procedure in the ET formulation can for example be found 
Ref. I BP 65|]. The rules for canonical LF-quantization have been taken from Refs. 



| CR+ 73i |CY 72 | 
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normal coordinates 



light-front 



coordinates 



X time 
x^,x'^,x^ space 



X 



+ _ x^ + 



^V2 



time 



X —- — -f^^x\x space 



scalar product 



a ■ b — 
= a^feo - ab 



a ■ b — 

a+b- + a-b+ - a^b^ - a^b^ 

— * 

= a'^b~ + a'b'^ — a±b± 



Lagrangian density 

c = i {docpf - 1 (v 0) ' - v{cj>) c = - 1 (^v± 0) ' - 1^(0) 

conjugate momenta 



canonical commutation relations 



[tt{x, t),ip{y,t)] 
— —iS^{x — y) 



[k{x ,x^,x+),ip{y ,|/±,.x+)] 
= -I5{x- -y-)5'^{x^-yA_) 



Hamilton operator 



= Jd^xn{(p,7r) 
H. — Trdof — C 



P+ = jdx' j£x±n{(p,7r) 
H — 7rd+(p — L 



= ^ (Px TT 



PVn>=^n|^n> 

P fixed 



momentum operator 

\/(p P_ = jdx~(Px± nd^cp 

P± — J dx~d^x±_ 7rd±(p 
eigenvalue equation 

P+\-^Jr,>= P+n|^n> 

P_,Pi fixed 



hadron masses 



M^ = E^- P2 

n n 



= 2P+„P_ - P2 



Table 2.1: canonical quantization in ordinary coordinates and on the light- 
front 
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dynamics and can start solving the equations of motion. Typically, one ei- 
ther makes some variational ansatz or a Fock space expansion. In the latter 
approach one writes the hadron wave function as a sum over components 
with a fixed number of elementary quanta (for example in QCD: gg, qqqq, 
qqg, e.t.c). The expansion coefficients, i.e. the wavefunctions for the cor- 
responding Fock space sector are used as variational parameters. They are 
determined by making the expectation value of the energy stationary with 
respect to variations in the wavefunction. Typically the variation is done for 
fixed momentum.0 This whole procedure results in coupled integral equations 
for the Fock space components. In general they have to be solved numeri- 
cally. In practical calculations, since one cannot include infinitely many Fock 
components, one has to introduce some ad hoc cutoff in the Fock space. Thus 
it is very important to demonstrate that physical observables do not depend 
on how many Fock components are included. 

Until one selects the canonically conjugate momenta and postulates equal 
r commutation relations, i.e. at the level of the classical Lagrangian, the 
transition from ET to the LF consists of a mere rewriting. After quantization, 
the independent degrees of freedom consist of the fields and their conjugate 
momenta on the initial surface = for ET and = for LF). Thus 
different degrees of freedom are employed to expand physical states in the 
ET and in the LF approach. Of course, after solving the equations of motion, 
physical observables must not depend on the choice of quantization plane. 
However, it may turn out that one approach is more efficient (e.g. faster 
numerical convergence) than the other or more elegant and more easy to 
interpret physically. In general, this will of course depend on the details 
of the interaction. An extreme example is QCD^^^(A''(7 — oo). In the ET 
approach [PG 78| , [Li 86| , |HN+ 88|] one first has to solve coupled, nonlinear 
integral equations with a singular kernel to obtain the Hartree-Fock solution 
for the vacuum. Then, in order to calculate meson masses, one has to solve 
the two body equation in this background, which amounts to solving another 
set of coupled (linear) integral equations with singular kernel. In the LF- 
approach [[bio 74| , |Ei 76|| all one has to do is solve one linear integral equation 
with singular kernel. The numerical results for the meson spectrum are in 
extremely good agreement between the two approaches, but numerically the 
LF calculation is more than one order of magnitude faster! In this case the 
simplification arises because the LF-vacuum is trivial — a point which will 
be elaborated in more detail below as well as in Chapter ^. 

Which approach is preferable may, however, also depend on the observ- 



•^On the LF this is very important because P+ oc 1/P-, i.e. unrestricted variation (P_ 
allowed to vary) results in P_ — > oo. 
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ables in which one is interested. The most prominent example is deep in- 
elastic scattering. As discussed in the introduction, parton distributions are 
much more easily accessible on the LF than in usual coordinates. 



2.2 e-Coordinates on Finite Light-Front 
Intervals 

One issue one may be worried about is the question of equivalence between 
the LF approach to field theories and other approaches. On the LF one 
imposes commutation relations at equal LF-time, i.e. between two space- 
time points that are connected by a light-like distance. Thus it is a priori not 
clear whether the initial value problem with initial conditions on a null plane 
is well defined ||Ro 70| , P^M 92|| and whether there arise any confiicts with 



causality on the LF. The situation becomes particularly worrisome when one 
introduces a "box" in the longitudinal x~ direction (to keep IR-singularities 
under control) and imposes periodic or quasiperiodic boundary conditions 
at the ends of the box — i.e. one imposes boundary conditions between 
points that may be causally related. One way to address this issue in a well 
defined way is to define the LF via a limiting procedure by starting from a 
spacelike quantization surface and carefully rotating this surface until one has 
'reached' the LF (note: although there are some similarities, this should not 



be confused with a Lorentz boost to infinite momentum |[We 69| , PR-|- 73| ). 

In order to be able to control infrared singularities, let us formulate the 
dynamics on a finite LF interval with extension L in the x~ direction. 
On a finite interval, boundary conditions have to be specified, e.g. 0(x~ + 
L,x~^) = (f){x~ ,x~^). However, if one is working on the LF, imposing boundary 
conditions means relating fields at points that are separated by a light-like 
distance — obviously one may run into trouble with causality at this point. 



To avoid this dilemma, Lenz et al. |[LT 91|| have introduced £-coordinates 
which are defined as follows, R 



x^ = X 



x+ = + (2.3) 



Now points at opposite ends of the interval (with coordinates {x^ +L, xf) and 
{x~,xf) are separated by a spacelike distance ds"^ = —2eL and no confiict 



^To simplify the notation, only 1+1 dimensional examples will be discussed in this 
section. 



See also Ref. pF Sgj . A shghtly different approach, where both x+ and x are rotated 



away from the Hght-cone, has been studied in Ref. |Ho 92 1 
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with causality arises from imposing boundary conditions. In e-coordinates 
the scalar product is given by 



A-B = A+B_ + A_B+ + 2-A+B^ 

1j 



(2.4) 



and thus the Lagrangian density (for the rest of this section, the subscript e 
will be dropped to simplify the notation) for </)f^^ reads 



(2.5) 



Since C is quadratic in (9+0, quantization in e-coordinates is straightforward 
(as in usual coordinates). One finds ||LT 91| 



n 



5C 



+ 

1j 



(2.6) 



n(a; ,x~^),(f){y ,x^) = —iS{x —y ) 



and 



H 



+ 



4! 



(2.7) 
(2.8) 



In these coordinates, the free dispersion relation (A = 0) is given by 



(2.9) 



where p-{n) = ^ as usual in a box with periodic boundary conditions. 
Later, we will also need the normal mode expansion of the fields with periodic 
boundary conditions 



[X 



E 



71" '^x/^n 



n(x-) = y: 



n 



(2.10) 



where u;^ = L^p_{nY + "^m? and the a, a) satisfy the usual commutation 
relations, e.g 

(2.11) 



The most significant difference between the dispersion relation in e- 
coordinates ( p. 91) and the dispersion relation on the LF (p+ = is the 
appearance of two solutions of p+ for each p_ in e-coordinates, while the 
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Figure 2.1: Free dispersion relation in e-coordinates versus the dispersion 
relation in the LF limit. 

dispersion relation on the LF yields just one solution for each p_ (Figure 
|2.1| ). For positive energy (p+ > 0) modes, the LF momentum p_ is positive 
whereas the momentum p_ in e-coordinates can be both positive and nega- 
tive. This has importance consequences for the vacuum structure which will 
be discussed in Chapter ^. 

In the limit —>■ {L fixed) the LF is recovered: 

n^9_0 (2.12) 

H I dx-'^^P' + (2.13) 

For all nonzero e, the relation between the momenta and the fields 
(|2.6|) contains the time derivative of the fields and the fields are quan- 
tized as usual (|2.7| ). However, for e = 0, Eq.( p.6|) becomes a con- 
straint equation, and the Dirac-Bergmann algorithm (see Appendix 0) yields 
[d-(f){x-, x+), (f){y-, x+)] = |5(x- - y-). 

It should be noted, that the order of limits does matter, i.e. it is im- 
portant whether one takes the LF limit (-^ — > 0) first or the continuum limit 
(L oo). This will be discussed in detail in Chapter |^. 

2.3 Examples for Canonical Light-Front 
Hamiltonians 
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2.3.1 Scalar Fields 

Self- interacting Scalar fields in the LF framework have been discussed in 
Refs. [|CR+ 73| , |CY 72|] . In order to keep the discussion as general as possi- 
ble, we will work in transverse dimensions, where Z)_l = 0,1,2. For a 
polynomial interaction ^, 



C= -da 
2 



m 



-jC 



int 



(2.14) 



where C is a polynomial in 0, the momenta conjugate to are 



n = 9_ 



(2.15) 



with commutation relations 



n(x ,fj_,x+),0(?/ ,y±,x^)\ = --S{x -y )5{x±-y±). (2.16) 
Note that this implies nonlocal commutation relations for the field 0, e.g. 



{x ,x±,x^),(j){y ,y±,x'^)\ = --e{x -y )6{x±-y±), (2.17) 

where e{x) = 1 for x > and e{x) = —1 for x < 0. The Hamiltonian 
density (P+ = / dx~ d^^xxTi) is obtained from Eq.( |2.14] ) via a Legendre 
transformation 



1 / ^ \ 2 



(2.18) 



The commutation relations ( |2.16| ) are easily satisfied if we make a mode 
expansion 



^fWkZJ (27r)^i/ 



(2.19) 



where ^ satisfy the usual boson commutation relations, e.g. 



6{k^ - qJ)6{k±_ - gl). 



(2.20) 



^In 3+1 dimensions, rcnormalizability restricts the interaction to 4th order polynomials, 
but in 2 + 1 or 1 + 1 dimensions higher order polynomials are conceivable (6th order and 
oo order respectively). 
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Longitudinal and transverse momentum operators contain no interaction 
terms 



dx~d^^x±UV_ 



dk. J d^'-kj^al^^^a.^^,^ (2.21) 



where normal ordering terms have been dropped. Most of the numerical 
studies of self-interacting scalar fields have been done in 1 + 1 dimensions 
[ [HV 871 |Bu 93|| using discrete light-cone quantization (Section ^3]) . A more 
recent work employs Monte Carlo techniques to solve c/f'^-theory in 2 + 1 



dimensions [Pu 94a| 



Complex scalar fields can always be reduced to real scalar fields by work- 
ing in a Cartesian basis ^ = (0i + ^02 ) and thus need not be discussed 
here. 



2.3.2 Fermions with Yukawa Interactions 

To keep the discussion as general as possible we assume an interaction of the 
form ipTipip, where is either scalar or pseudoscalar and F is either 1 or 275 

£ = ^(i^-M - gT(j)) + ^ (df,(pd^(p - m^^^) . (2.22) 

One novel feature compared to normal coordinates and compared to self- 
interacting scalar fields on the LF is the fact that not all components of i/j 
are independent dynamical degrees of freedom. To see this, let us introduce 
projection matrices 'p(='=) = ^7^7^ where 7^ = (7° ± 7^) / ^/2. Note that 
7+7"*" = 7~7~ = implies = p(-)'p(+) = 0. These projection 

matrices can be used to decompose the fermion spinors into dynamical and 
non-djTiamical components = ^/'(+) -|- ?/'(_) , where '?/'(±) = V^^^il). The 
Lagrangian does not contain a LF-time derivative (9+) of ■?/'(_) 

C = 72^1^)29+ + y2V'(-)«5-^(-) - + 7°-M) ^(-) 

-^f_) [laj^ + 7°-M) ^(+) + \ (5^05^0 - m^'') , (2.23) 

where d± = 7'^7_l and A4{x) = M + gr(f){x). Thus the Euler-Lagrange 
equation for ?/'(_) is a constraint equation 

y2i9_^/'(_) = {ia±d± + 7°A^) ^(+). (2.24) 
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It is therefore necessary to eliminate the dependent degrees of freedom 
before quantizing the theory. Here we proceed by solving Eq.( p.24 ) and 
inserting the solution back in the Lagrangian ( p.23| ), yielding [ pY 72|| [] 



C 



(+) 



-— {taJi_ + 7°A^) ^ {tdJi_ + 7°-M) ^(+). (2.25) 



V2 



The ambiguities associated with the inversion of the differential operator will 
be discussed in Section p.2.2| . Here we just define 

dy~e{x'' - y')fiy~,x±). 



(2.26) 



The rest of the quantization procedure is now straightforward. The Hamil- 
tonian is given by 

P~ = j dx-d^^x±n (2.27) 



where 



V2 

1 

2 



+ m 



2i2 



(2.28) 



Note that Eq. ( |2.28| ) contains four-point interactions of the form 
(^t^-)~^ 4>i^i+) which were not present in the original Lagrangian (|2.22|) . 
Note also, that the fermion mass M enters the Hamiltonian density ( ^.2^ ) 
in two different places: in the kinetic term for the dynamical fermion 
field, Hkin oc M^?/'|^_(_^ (z(9_)~"^ as well as in the three point vertex, 

Mgipj^^^ (^f^-)~^ r0V'{+) + h.c. In Chapter ^ we will find that in general 
these two masses are renormalized differently. 

The scalar field is quantized as in Section ( p.3.1| ). For the fermions, one 
imposes anti-commutation relations only for the independent component 



(2.29) 



with {V^(+)(a;),^(+)(l/)}^^^^^ and {'^{+){x),ip{+){y)} ^^^^^ both vanishing. 

^Another option would be to use Dirac-Bergmann quantization (Appendix |A|). Up to 
possible differences in the zero-mode sector, the result is the same. 
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For practical calculations it is very useful to make a mode expansion. 
Let u{p,s), v{p,s) be the usual particle and antiparticle spinors, satisfying 
{'^—M)u{p,s) = and {-^ + M) v{p, s) = 0, where s labels the spin. The 
normalization is fixed such that 

a/2m{_^)(p, s) = u{p, s')7-m(p, s) = 2p_5ss' 

v^^(+)(p>s>(+)(p, s) = v{p,s')-f-v{p,s) = 2p^6ss'- (2.30) 

For ■?/'(+) we make a plane wave ansatz 

^(+)(^) = / (27r)^x'/2 I E [bp-P^sUi+){p, s)e-'P'' + 4_p.^,i;(+)(p, s)e 

(2.31) 

One can easily verify that in Eq.( |2.3lD satisfies the anti-commutation 
relations above ( p.29| ), provided 

{bl_p^r^ bq-<i±s} = 5(p- - q-)S{p± - q±)6rs 

{4_Pxr'Vg±4 " ^{P- - (1-)^{P± - Q±)^rs (2.32) 

with all other anti-commutators vanishing. 

Nonperturbative numerical works on the LF Hamiltonian with Yukawa 
interactions (|2.28|) have been restricted to DLCQ calculations in 1+1 di- 
mensions |[PB 85|| as well as to 1+1 dimensional |pP 91|| 3+1 dimensional 
[ |Wo 93|| calculations which use Tamm-Dancoff truncations to fermion and at 
most two bosons or fermion, antifermion and at most one boson. 



ipx 



2.3.3 QED and QCD 

Before one can canonically quantize a gauge field theory, one must fix the 
gauge — otherwise one has to deal with the infinite degeneracy associated 
with the gauge symmetry. 

In the context of LF quantization one usually picks the LF-gauge A'^ = 
y4_ = 0. There are several (related) reasons for this choice. In QED, the 
constraint equation for the "bad" spinor component reads 

v^iD_^(_) = V2{id^ - eA_) = {ia±D^ + 7°m) (2.33) 

The solution to this constraint equation, 

^(_) = ^ {id. - eA.)-' (ta^D^ + 7°m) ^(+), (2.34) 
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contains in the denominator and, unless one chooses = gauge, one 
thus obtains terms which have A_ in the denominator in the LF-Hamihonian. 
In other words, in any gauge other than the LF-gauge the canonical LF- 
Hamiltonian always contains all powers of A_ (after expanding the geometric 
series) appearing in the interactions. 

In QCD one faces the additional problem that II^, the momentum conju- 
gate to A_i_, satisfies a nonlinear constraint equation if A_ ^ (Section |3.2.3| ). 
Another reason to pick LF gauge is that A_ = is invariant under the kine- 
matic Lorentz symmetries of the LF, i.e. under all transformations that leave 
the plane = invariant. It is for these reasons that the LF gauge has been 
commonly used for canonical LF quantization of gauge field theories, and will 
also be used here — despite all the difficulties which are inherent to the LF 
and axial gauges ||Ma 83| , |BD+ 85| , |Wu 77| , |LN+ 94a| , |LN+ 94b| , |LN+ 94c| . 



Even after fixing the gauge, not all degrees of freedom are dynamical 
(similar to their time derivative does not enter the Lagrangian). Before 
we can proceed with the canonical quantization we first have to eliminate 
these dependent variables by solving those equations of motion which are 
constraint equations. For we use Eq.( |2.34|) (note: A_ = 0) and proceed 
similar to the example of the Yukawa theory. Since the time derivative of A+ 
does not enter the Lagrangian, has to be eliminated as well, by solving 
its constraint equation (obtained by varying the Lagrangian density with 
respect to A^) 

d^_A+ = d-V±A^-j+, (2.35) 
where j+ = v^e?/'|_^ -,'?/'(+), in QED and 

= 9_ Vx^a + gr'^A^bd-A^c - J+, (2.36) 

where = v^fi'^/'(_,_)Q,V'(+)/3^! iii QCD. After inserting Aj^ back into the 
Lagrangian one can proceed with the quantization as usual. One finds 

1 1 
Hqed = ^^f+) + 7°m) — (za^Dx + 7°m) ^(+) 



- [d-V^A^ - j+ ) ^ (9_ V^lx - j+) , (2.37) 



and 



Hqcd = -^^f+) + 7°m) ^ (m^D^ + 7°m) 

-I {D^d^A^a - ^ {D^d^A^, - j+) (2.38) 
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where D±d-A±a = ^±d-A±a + gf"''"^A±bd-A±c- The commutation relations 
are similar to the ones in Yukawa theory 



QED: 



V2 



[d.AUx),A^j{y)l^^ 

QCD: 



:y+ = -2^(^ - y )^(^^ - y^)^ij 



V2 



V^'^^6{x -y )S{x± - y±)6ai3 



[d-A±aiix),A±bjiy)] 



-5{x - y )S{x± - y±)6ij6a 



(2.39) 



Similar to the approach to scalar field theories and Yukawa theories, one 
may now attempt to solve the above Hamiltonians by making a mode expan- 
sion and using matrix diagonalization. In 1 + 1 dimension this method was 
very successful ||EP+ 87| , [EF^ |HP+ 9^ , |Bu 89ai |Bu 89b| , [BB~gT| . In 3 + 1 



dimensions, this approach suffers from a fundamental problem 0: charged 
particles are subject to a linear, confining interaction — which is present 
even in Hqed- For gauge invariant amplitudes (all intermediate states in- 
cluded, which contribute to a given order of the coupling) this linear potential 
is canceled by infrared singular couplings of charges to the ^-components of 
the gauge field. However, in most practical calculations, drastic truncations 
of the Fock space are used to keep the dimension of the Hamiltonian ma- 
trix within practical limits [[TB 91| , |KP 92|| . This approximation results in 
incomplete cancelations of IR singularities and IR divergences result. Partly 
responsible for this disaster is an improper treatment of zero-modes and in- 
complete gauge fixing. If one integrates the Maxwell equation for F^^ over 
x~ one finds [[KP 93|| 



— dx y4_(x^,a; ,x±) = dx j^{x^^x ,x±) 



(2.40) 



i.e. in general, when / dx~j^{x^,x~,x±) ^ 0, the "gauge" A_ = is in- 
consistent with the equations of motion. On a finite interval, with peri- 
odic boundary condition, this becomes clearer because then a Wilson loop 
"around the torus", exp (ie_f dx~ A^{x^ ,x~ ,xsy), is a gauge invariant quan- 
tity. The closest one can get to the LP gauge is d-A_ = 0. In this gauge 
one can now investigate the problem of incomplete gauge fixing. The gauge 
d^A^ = still leaves the freedom of ^"-independent gauge transformations 
A' = A^ + d^x where (9^x = (or = if we restrict ourselves 



Besides numerical difficulties which will be discussed in Section 5.1 
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to periodic x) ||KP 93|| . In such an incompletely gauge fixed situations, not 
all degrees of freedom are physical and approximations may result in incon- 
sistencies. A typical example is the residual or transverse GauB' law ( |2.4CI| ), 
which is a constraint on the physical Hilbert space. Such constraints must 
either be imposed on the states or one can also use them to eliminate "un- 
physical" degrees of freedom (here dj_ J dx~A^). The abovementioned, in- 
complete cancelation of IR singularities in the Tamm-Dancoff approximation 
occurs because the transverse Gaufi' law ( |2.40|) is violated. A more thor- 
ough discussion on this subject and possible caveats can be found in Refs. 
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Chapter 3 

The Light-Front Vacuum 



3.1 The Physical Picture 

In the Fock space expansion one starts from the vacuum as the ground state 
and constructs physical hadrons by successive apphcation of creation op- 
erators. In an interacting theory the vacuum is in general an extremely 
complicated state and not known a priori. Thus, in general, a Fock space 
expansion is not practical because one does not know the physical vacuum 
(i.e. the ground state of the Hamiltonian) . In normal coordinates, particu- 
larly in the Hamiltonian formulation, this is a serious obstacle for numerical 



calculations. As is illustrated in Table the LF formulation provides a dra- 



matic simplification at this point. While all components of the momentum in 
normal coordinates can be positive as well as negative, the longitudinal LF 
momentum P_ is always positive. In free field theory (in normal coordinates 
as well as on the LF) the vacuum is the state which is annihilated by all 
annihilation operators a^. In general, in an interacting theory, excited states 
(excited with respect to the free Hamiltonian) mix with the trivial vacuum 
(i.e. the free field theory vacuum) state resulting in a complicated physical 
vacuum. Of course, there are certain selection rules and only states with the 
same quantum numbers as the trivial vacuum can mix with this state; for 
example, states with the same momentum as the free vacuum (P = in nor- 
mal coordinates, P_ = 0, Pj_ = on the LF). In normal coordinates this has 
no deep consequences because there are many excited states which have zero 
momentum. On the LF the situation is completely different. Except for pure 
zero-mode excitations, i.e. states where only the zero-mode (the mode with 
= 0) is excited, all excited states have positive longitudinal momentum 
P„. Thus only these pure zero-mode excitations can mix with the trivial LF 
vacuum. Thus with the exception of the zero-modes the physical LF vacuum 
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normal coordinates 



light-front 




t ^ _ m2+fe2 



k-,k± 

vacuum (free theory) 



a,|0) = 



ait_,fex|0) = 



vacuum (interacting theory) 



many states with P = 
(e. g. aialj^l 



|0> very complex 



fc_ > 
only pure zero-mode 
excitations have P_ = 

|0> can only contain 
zero-mode excitations 



Table 3.1: Zero Modes and the Vacuum 
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(i.e. the ground state) of an interacting field theory must be trivial.^ 

Of course, this cannot mean that the vacuum is entirely trivial. Otherwise 
it seems impossible to describe many interesting problems which are related 
to spontaneous symmetry breaking within the LF formalism. For example 
one knows that chiral symmetry is spontaneously broken in QCD and that 
this is responsible for the relatively small mass of the pions — which play 
an important role in strong interaction phenomena at low energies. What it 
means is that one has reduced the problem of finding the LF vacuum to the 
problem of understanding the dynamics of these zero-modes. 

First this sounds just like merely shifting the problem about the structure 
of the vacuum from nonzero-modes to zero-modes. However, as the free 
dispersion relation on the LF, 

- (3.1) 

indicates, zero-modes are high energy modes! Hence it should, at least in 
principle, be possible to eliminate these zero-modes systematically giving 
rise to an effective LF field theory [p^T 91 ]. 



Before we embark on theoretically analyzing zero-modes, it should be em- 
phasized that zero-modes may have experimentally measurable implications. 
This is discussed in Refs. ||Bu 92c|, [Bu 95| . 



3.2 Examples for Zero Modes 

Usually, in the context of LF quantization, fields that do not depend on 
x~ are called zero- modes (regardless whether they depend on x_l or not). 
However, for practical purposes, the following classification scheme seems to 
be particularly useful [[KP 93|| : If one denotes 

{f)o = ^f^dx-f{x-,x^), (3.2) 

then 

(/) - (^/^'^^(/)«= ^Z(^/^'^^X/^"/(^~'^^) (3-3) 
is called the global zero-mode, while 

h {f)o - (/) (3.4) 



^ Cases where the LF Hamiltonian has no ground state will be discussed below. 
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is called proper zero-mode. The "rest" , i.e. 

n 

/=/-(/)o (3.5) 

is called the normal mode part of /. The motivation for this distinction arises 
primarily from the fact that usually only the global zero-mode can develop a 
vacuum expectation value but also since proper and global zero-modes have 
a very different dynamics. 

Zero modes occur in various contexts and it is not yet entirely clear to 
what extend the various zero-mode effects, which will be discussed below, 
are connected. 



3.2.1 Constant Scalar Fields 

In (f)^ theory. 



C = -d^cl^dy - —<f - (3.6) 



if one chooses the "wrong" sign for the mass {m? < 0, A > 0), spontaneous 
symmetry breaking occurs already at the classical level. The field develops 
a vacuum expectation value and the symmetry — > —0 is spontaneously 
broken. At least in 1 + 1 and 2 + 1 dimensions, with appropriate values for 
the renormalized mass, a similar behavior is observed in the quantum version. 
Clearly, such a scenario requires a zero-mode. In the case of 0^ theory, one 
may imagine that a redefinition 

0-^0+(O|0|O) (3.7) 



eliminates the VEV of the global zero-mode [[HV 87|] . However, this does not 
mean that one has eliminated the zero-modes. In fact, by integrating the 
equations of motion over x~, one finds 

O = m2(0) + A±|'^rf^-03^ (3.8) 

which relates the zero-mode part of the field to the normal mode part. 
Clearly, this nonlinear operator identity implies that (for finite L), a mere 
shift of the scalar field is not sufficient to completely eliminate the zero- 
mode. Instead, two main classes of approaches are being used get the 
zero-modes under control. In DLCQ one attempts to solve the zero-mode 
constraint equation [Eq. ( |3.8| )] using various approximation or expansion 
schemes |PV 93| , PV 94| , PV+ 95|] . Due to nonlinear effects and opera- 



tor ordering ambiguities, solving Eq.(^.8[) becomes a nontrivial endeavor 
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PV 93| , PV 94 PV+ 95| , pF93| , PK+ 92a|| . In the other approach (the ef- 



fective LF-Hamiltonian approach, which will be discussed in detail in Sec- 
tion 3.3) one makes use of the fact that zero-modes freeze out for L 00. 
Instead of keeping zero-modes explicitly, one allows for an effective Hamil- 
tonian, which should account for their effects on normal modes in the large 
volume limit Pu 93|| . 



So far it is not known whether either one of these approaches to LF- 
quantization (explicit zero-modes and effective LF-Hamiltonian) leads to a 
consistent formulation of (p^ theory in the broken phase. It is also not known 
to what extend the particle spectrum in the equal time formulation agrees 
with the spectrum on the LF. Since the broken phase of 0^ in 1-1-1 dimensions 
has a rather rich spectrum: mesons, solitons 0, bound states and scattering 
states in the soliton-antisoliton sector 0, this seems to be an ideal test case 
for the various approaches to scalar zero-modes on the LF. So far, all works 
on (j)'^ on the LF have concentrated on demonstrating that spontaneous sym- 
metry breaking occurs and on reproducing the numerical value of the critical 
coupling constant from ET quantization.^ 

One of the most striking consequences of spontaneous symmetry breaking 
in (pfj^i is the emergence of solitons. While most LF workers choose bound- 
ary conditions that make it impossible to study solitons, soliton-antisoliton 
scattering states are often still possible. These states often have a very clear 



signature ||Bu 93|| and one can easily determine their threshold. Considering 



the extensive literature on LF-0f^^ (see e.g. Refs. | |FV 93| , [FV 94| , ^\/+ 95 



and references therein), it is surprising that solitons have been ignored so far. 

3.2.2 Fermionic Zero Modes 

Consider the free Dirac equation 

= (i7^9^ - M) ^ = (i7_9+ + i7+9_ - i^±d± - M) ip. (3.9) 
Multiplying Eq. (|3.9|) with V^'^\ where p(='=) = 7^7=*= /2 are the projection 



matrices introduced in Section 2.3.2, one obtains 



z7-(9_^/'(") = {i^±d± + M) (3.10) 

with i)^^^ = V^^H. Clearly Eq. (|3l0|) is a constraint equation and one must 
eliminate before one can canonically quantize the theory (the kinetic 



^Often excluded by boundary conditions on the fields. 
■^In general not excluded by boundary conditions. 



^In view of the nontrivial renormalization effects on the LF (see Chapter Q), comparing 
critical coupling constants on the LF with those from ET quantization is very treacherous. 
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term in the fermionic Lagrangian does not contain a LF-time derivative of 
ip^~^)- For all modes but the zero-modes this is straightforward. However, 
Eq.( |3.10| ) does not determine the ^"-independent components of ip^~\ In 
other words, because of possible "integration constants" , there is some am- 
biguity in defining the inverse of the differential operator 9_. 

For scalar fields, the time derivative is always accompanied by a space 
derivative (kinetic term: (pdj^d^cj)). Therefore, the zero-mode for scalar fields 
is not a dynamical degree of freedom, since its time derivative does not enter 
the Lagrangian. For Dirac fields this is different, since there the Lagrangian 
is linear in the derivatives, and the fermionic zero-mode is a dynamical degree 
of freedom. Little is known in this case beyond perturbation theory (see e.g. 

Refs. prmimgi). 



3.2.3 Gauge Field Zero-Modes 



For practical reasons one would like to work in the LF gauge = when 
quantizing gauge fields on the LF. The reason is that, only in the LF gauge 
are canonical field momenta simple. For example, in QCD, the kinetic 
term for the gauge field in the Lagrangian, —^F^uF^^ contains terms like 

-D+, Ax D_,A± , i.e. in general, the term multiplying d+A± contains in- 
teractions. As usual in LF coordinates, the canonically conjugate momentum 
satisfies a constraint equation 



n, 



6d^A, 



d-A± - ig 



(3.11) 



Only in the LF gauge is the constraint equation for linear in the fields, 
and one obtains simple commutation relations between the fields. 

The problem with the LF gauge, as with axial gauges in general, has to 
do with infrared singularities, particularly in the nonabelian case. In order 
to arrive at a well defined formulation of the theory, it is often very helpful 
to formulate the theory in a finite 'box' with periodic boundary conditions 
(i.e. a torus). That way, it is generally easier to keep track of surface terms 
that appear in formal manipulations which include integrations by parts. 

If one starts from an arbitrary gauge field configuration on a torus, it 
is in general not possible to reach the LF gauge (or spatial axial gauges) 
by means of a gauge transformation | Ma 85 , KP 93|| . This can be easily 
shown by considering the Wilson loop around the torus in the x~ direction: 
W = Pexp{ig § dx~A_). This is a gauge invariant quantity and thus does 
not change under a gauge transformation. If it were possible to reach the 
LF gauge, A^ = 0, be means of a gauge transformation this would mean 
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transforming W to 1, which is a contradiction. It turns out that on a torus, 
the closest one can get to the LF gauge is d_A_ = 0, i.e. the zero-modes for 
A_ remain and, due to their relation to the Wilson loop around the torus. 



they have a gauge invariant meaning [[KP 93|| . They are dynamical degrees 



of freedom (their derivative enters the Lagrangian). The zero- modes of a* 
behave very similar to a scalar field, in the sense that their time derivative 
does not enter the Lagrangian and hence they are not dynamical degrees of 
freedom. Recently, Kalloniatis, Robertson and collaborators ||KP 93|| have 
developed a systematic scheme to disentangle and resolve the various zero- 
mode problems that appear in QED and QCD. For example, projecting the 
QED Maxwell equations onto the proper zero-mode sector, they obtain: 

-dlA^ = gf (3.12) 
-2dl a' -dl a' -2d,d+ A = gdl T (3.13) 
-dl A +did+ A^ +d,dj A = gdl J , (3.14) 

where is the fermionic current operator. The first of these equations ( p.l2| ) 
is a constraint equation and can be used to eliminate the proper zero-mode 
of A'^ in terms of the current Q 

o+ 1 + 

A = -9-^ J , (3.15) 

which again demonstrates that A'^ = is in general not consistent with the 
equations of motion. 

Further simplification can be obtained by taking the (transverse) diver- 
gence of Eq. (|3.14| ) , yielding 



d+A'=g^d,J . (3.16) 



Inserting this back into Eq.( 3.14j) , one finds 




^l{5^,-^] A=g{Si,-^] J , (3.17) 



which can be used to eliminate the transverse projection of the proper zero- 
mode of A^ . Note that so far we have not yet completely fixed the gauge, 
since d^A^ sill leaves the freedom of purely transverse gauge transformations. 



^In the charge neutral sector, the global zero-mode of J"*" vanishes and thus the inverse 
Laplace is well defined. 
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A^^ — > + d'^Q, where fl = Q{x^,x±). One can use this residual gauge 

o « 

freedom to set di A = 0. In combination with Eq.( p.l7l ), this completely 



determines the proper zero-mode of A^. Up to this point, it seems that the 
zero-modes in QED pose no real problems in the LF formulation. 

The real problems in this formalism arise when one tries to implement 
these results in a quantum formulation that includes fermions. This can be 
seen when one inserts the solution for / dx~A_ back into the Lagrangian, 
yielding a four Fermi interaction term of the form 



^ ^ ■ <3-^**' 



Similarly, inserting the solution for A into the Lagrangian yields 

r _ djdj 

The presence of such terms, which contain the "bad" current j+ = 
\/2e'?/'(_-)V'(-) leads to nonlinear constraint equations for '?/'(-)• Because of 
the difficulties in solving this nonlinear constraint equation, it has so far not 
been possible to write down the LF Hamiltonian for QED or QCD in terms 
of physical degrees of freedom and including all zero-modes, in closed form. 
Only perturbative expressions for the Hamiltonian in terms of physical de- 
grees of freedom have been found so far ||KP 93|| . Similar problems arise in 



the DLCQ formulation of QCD with additional complications arising from 
the difficulties in quantizing the gauge field when A"*" ^ 0., arising from 
the nonlinear constraint relation between fields and their canonical momenta 



From the practitioner's point-of-view, it would be helpful to know to what 
extend this elaborate machinery is actually necessary if one is interested only 
in the large volume limit. On a finite interval, gauge field zero-modes clearly 
play an important role. For example, they are essential to generate the 
correct potential for a heavy quark- ant iquark pair in 1+1 dimensions on a 
circle in Coulomb gauge | [ES 95| , [KP+ 94| ]. However, in the latter example, 
zero-mode effects for color singlet states disappear in the limit of a large 
interval. Unfortunately, it is not clear whether this result carries through to 
higher dimensional gauge theories. 

3.2.4 Perturbative Zero-Modes 

The zero-modes discussed are either connected to purely nonperturbative ef- 
fects (like in the case of spontaneous symmetry breaking for scalar fields) 
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or seem to be at least connected with nonperturbative physics (hke infrared 
singular long range effects for gauge fields). There are, however, plenty of 
examples where zero-mode effects appear already on the level of perturbation 
theory. Examples include disconnected vacuum diagrams ||CM 69|| , "gener- 
alized tadpoles" for self-interacting scalar fields |Gr 92a| , |Bu ! 
"rainbow diagrams" for the fermion self-energy [ PL 91a| , |BL 91b|| 



as well as 
These 



examples will be discussed in more detail in Chapter |J. 

In perturbation theory in LF gauge the gauge field propagator 



k-n 



+ is 



(3.20) 



{n-A = AS) becomes singular as /c_ — 0. There exist several "prescriptions" 
to handle this singularity. The most useful prescription for perturbative 
calculations is the Mandelstam-Leibbrandt (ML) prescription ||Ma 83|| , where 
one replaces 



1 ^ ^ ML 1 
k-n /c_ k- + resign (A;_| 



k. 



k.k^ + ie 



(3.21) 



The crucial property of this prescription is that the pole structure is similar 
to the one of a typical Feynman propagator, with poles in the second and 
fourth quadrant of the complex fco-plane, and thus allows to perform a Wick 
rotation. This is not the case for the principal value (PV) prescription 



k. 




1 



k_ + ie 



IE 



(3.22) 



with poles in the first and fourth quadrant. 

One of the major disadvantages of the ML prescription is the fact that 
it introduces additional energy (A;+) dependencies in the propagator, which 
cannot be generated by a canonical LF Hamiltonian |[RM 94|| . However, re- 
cently the ML prescription has been successfully implemented in a LF Bethe- 
Salpeter approach to bound states |[LS 93|| . Conversely, in (^CY)y+i^^c oo) 
the ML prescription ||Wu 77| yielded a spectrum that disagreed with the 
canonical LF approach |Po 74|| as well as with the result from equal time 
quantization |P!ji 86 . More recently, light-like Wilson loops in 1-1-1 dimen- 
sions have been calculated, using various prescriptions for gauge field propa- 
gator ||BD-|- 94|] , and it was found that only the principal value prescription 
yields the exact area law one expects for gauge fields in 1-1-1 dimensions (on 
noncompact manifolds). 
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3.3 Zero Modes and the Vacuum in 
^-Coordinates 



3.3.1 General Considerations 



For a free particle P+ = ^ + \/p- + 2em? / l}j and p_ is no longer 
restricted to positive values (Fig. Therefore, for all finite values of e/ 



the vacuum in e-coordinates is nontrivial. Since e-coordinates (see Section 
|2.2| ) provide a controlled and well defined approach to the LF, it seems very 
natural to employ this framework for studying the LF vacuum. 

Let us first consider the canonical LF limit (L fixed, e ^ 0). In this 
case it is straightforward to derive an effective LF-Hamiltonian from the e- 
Hamiltonian [|LT 91|| (for a related work, see Refs. |PN+ 95| , [VF 94|| ). For 



finite L the momenta are discrete. Without interactions the energy of the 



zero-mode (p+(0) = my and the energy of modes with negative momenta 
(p_i_(— ra) ^ diverge as e — 0, while the energy of all positive momentum 
modes {p-{n) ^ remains finite. For interacting fields there will be 

some slight quantitative changes, but the general picture should remain the 
same: zero-modes and negative momentum modes are high energy modes 
— separated from positive momentum modes by an energy gap of O (y^) 

and O respectively. Thus although p_ < modes may acquire nontrivial 
occupations, p_ > modes have too little energy to cause any excitations 
within the p_ < sector for e — >■ 0: the p- < modes freeze out and can be 
replaced by their vacuum expectation value (VEV). 

At this point it seems that we have succeeded in deriving a nontrivial 
effective LF-Hamiltonian. Unfortunately, we arrived at this result by ap- 
proaching the LF in such a way that the invariant length of the interval 
(cx: Le) approaches zero, i.e. as discussed in Ref. ||LT 91|| , the effective theory 



that we have obtained is not necessarily equivalent to the original covariant 
theory. This can be easily illustrated by means of a perturbative example. 
Consider a simple tadpole with a mass insertion (to make it convergent) in 
1 + 1 dimensions 

cPk 1 i , . 

(3.23) 



(27r)2 (p -^2 _^^0)' 47rm2 



On a finite interval (with e coordinates) one obtains instead (/c_ (n) = fn 
^ 1 ^ r dk4 



L 27r (2ej^2^ ^ 2k+k^ - m2 + iO 
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-1 -3/2 



+ m 



(3.24) 



Clearly, in order to recover the continuum result ( p.23| ) one must take limits 
in such a way that the invariant length if the interval becomes infinite. If one 
takes the LF limit first (e — ^ 0, L fixed), one obtains a divergent contribution 
from the zero-mode. 

A different result is obtained if one performs the continuum limit first 
{L oo, e/L fixed). Since this corresponds to sL oo no problems with 
perturbation theory arises. However, since the spectrum is now continu- 
ous, there is no mass gap and the derivation of the effective Hamiltonian 
for e/L —>■ becomes more complicated. Nevertheless, it is still possible: 
first, note that the momentum scale of the continuum Hamiltonian is m^e/L 
since momentum dependent terms in the continuum Hamiltonian appear 
only in the kinetic term oc ^ — A;„ + y/c?. + m'^2e/L and in vertex factors 

oc + m?2e/ L^j. Thus the typical momentum scale in the vacuum is given 

by p'"_^^ = O (^^J^rnj . Similarly the energy scale for vacuum excitations (zero 



total P_) is of the order O y\J^rrvj. Suppose one is interested in the effective 
Hamiltonian for a physical particle of total momentum p^°^ moving in the 
vacuum. If |^ ^ 1 then there is almost no overlap between the wave function 



of the vacuum p"^"^ = O (y^^) and the wave function of the partons in 

the particle p^'*'"*"" because the parton wavefunction (calculated for example 
with a typical LF Hamiltonian) vanishes for small momenta. Thus one can 
introduce an energy gap by hand without affecting the dynamics in the limit 
■| — > 0: for example, by selecting cutoffs Ai and A2 such that 



m 



< Ai < A2 < 



(3.25) 



and removing all modes with Ai < fc_ < A2. First, this gives rise to a mass 
gap and one can argue that the modes with fc_ < Ai remain frozen (energy 
scale A;+ > ||-) when excitations with > A2 are present (energy scale 

k+ < 1^): in second order perturbation theory, the energy shift for modes 
with k_ > A2 due to excitations of n modes with k_ < Ai is given by 



oc 



Ai 



n/2 



-A^ = -Ai 



(3.26) 



Here A^ "''^ is a vertex factor, arising from the factor in the expansion of 
the fields (p (Eg. ( |2.1d| )), the factor A" is the phase space factor for n modes 
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with = 0{Ai) and states with < Ai excitations are off-shell by at least 
j_ 

Ar 

Since A;_ < Ai excitations are suppressed, the effective LF-Hamiltonian 
for the modes with A;_ > A2 contains the A;_ < Ai modes only via their VEV 
(which may be nontrivial!) 

V,'!U = {0k^<M\V\0k_<A,). (3.27) 

The crucial point is that the parton distribution calculated with such an 
effective LF Hamiltonian vanishes for small momenta in the above super- 
renormalizable example Thus as long as A2 is small enough compared to 
the total momentum of the particle p, the parton distribution vanishes al- 
ready for momenta much larger than A2 and the presence of the cutoff does 
not affect the parton dynamics. Since the VEVs are nearly independent of 
J, so is the effective Hamiltonian. Thus the suppression of the parton dis- 
tribution due to the kinetic energy sets in at a value cp^*, where c is nearly 
independent of j^. Thus, for |; — >^ 0, A2 can easily be chosen smaller than 
^ptot t^i^Wq ( p.25[ ) remains satisfied. In other words, A2 can be chosen such 
that the parton dynamics is independent of the exact position of the cutoff. 
Similarly, since vacuum momenta are restricted to p'"_^'^ = O (^^J^rrij <^ Ai, 
the presence of the cutoff does not affect the dynamics of the vacuum either, 
i.e., the numerical value of the VEVs which enter Eq.( |3.27| ) is independent 
of the cutoff. 

In the 2nd order perturbation theory argument above we made use of 
Ai ^ A2 to make sure that the energy denominator in Eg. ( p. 261) is of 



the order 0(^-^. This is actually not necessary, since the occupation of 
these modes is anyway dynamically suppressed for k_ < cp^ and as long as 
Ai -C cp_, the energy denominator will automatically be O (^j^^ or smaller. 

Thus we can actually let A2 — > Ai, i.e. remove the cutoff, without altering 
the conclusion. Introducing a mass gap was helpful in deriving an effective 
Hamiltonian for modes with j9_ ^ "^y^- However, since the solutions of the 
effective Hamiltonian vanish at small p- anyway, there is no need for a cutoff: 
a region void of excitations between rn^J^ and cp- develops dynamically 
(Figure p.l|) and this is sufficient to derive an effective Hamiltonian. 

In the end, the following result is obtained. Suppose we started from 
some polynomial interaction 

£-* = -^A„^. (3.28) 



^Roughly speaking, the LF kinetic energy T, which one can calculate from the parton 
momentmxi distribution f(k^), using T = J^- dk^ , has to remain finite. 
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Figure 3.1: Schematic occupation of modes in the presence of a particle 
with momentum p_ for e/L -C 1 (i.e. "close to the LF"). The modes near 
k- = are already present in the vacuum and are dynamically restricted 
to = O {m{e / L)'^^/'^^^ The "parton distribution", i.e. the modes which 
are occupied in the presence of the particle but not in the vacuum, vanish 
at small k- at a momentum scale which remains finite as e/L 0. The 
presence of the cutoffs has almost no effect on the dynamics. 

Then, using Eg. (|3.27|) (after some combinatorics) the effective interaction, 
which enters the LF Hamiltonian for p_ > modes in the limit |; ^ is 
given by 

^r// = - E An E 7^^(01^10) (3.29) 

(in order to obtain this result one also uses that, after normal ordering, the 
p_ > modes do not contribute to the VEVs). Eq.( |3.29| ) is a remarkable 
result. It states that nontrivial vacuum effects enter the LF-Hamiltonian 
only via effective interactions. The effective coupling constants depend on 
the vacuum condensates which, in general, cannot be obtained directly from 
a LF calculation Q They must be considered as renormalization parameters 
of the LF theory. Eg. ( |3.29D is also valid in situations where spontaneous 
symmetry breaking occurs. For example in 0^ theory, (O|0|O) may become 
nonzero and a 0^ interaction will thus appear in the effective Lagrangian. 
However, note that only a finite number of condensates is necessary to 
specify the effective LF Hamiltonian: if is the highest power of entering 

^However, there are exceptions where one can use sum rules or consistency conditions 
to determine the effective coupHngs iteratively. Examples will be discussed in the following 
section. 

^ It should be emphasized that we did not make any mean field assumptions, such as 
(0|(?!)'^|0) — (01010)'', in order to arrive at this result. 
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the canonical LF Hamiltonian then only condensates (010*^10) with k < N 
need to be considered. 

At several points in the above discussion it was important that the theory 
is free of divergences (up to a finite number of diagrams which can always be 
subtracted before applying above argumentation). First this was important 
to insure that the momentum scale in the vacuum is finite. Secondly, it 
was important because only in the absence of divergences one can apply the 
kinetic energy argument to prove that the parton momentum distribution 
vanishes for small parton momenta. Therefore one must be very careful when 
generalizing the above results to higher dimensional field theories. Eq. ( |3.29| ) 
can be used in renormalizable field theories only after a cutoff has been 
imposed. This is for example the case for the transverse lattice which will be 
discussed in Section pl3| . 



3.3.2 A simple Example for the limit L 



oo. 



e/L^O 



The appearance of the gap as L — > cxd (first) and e/L — > is best understood 
by studying a concrete example. Ideally this implies considering some non- 
trivial interacting field theory and calculating the occupation of the modes 
nonperturbatively for various e/L and L —>■ oo. However, even in integrable 
models, such as the sine-Gordon model, the occupation of the modes is not 
known exactly! Since numerical calculations at small but finite e/L and 
L — > cxD are very complicated — particularly if one is interested in momenta 
of the order of m^e/L — and we will proceed by studying a perturbative 
example. Due to the fact that the appearance of the gap is mostly a conse- 
quence of dimensional analysis, this will be sufficient to highlight the essential 
physics of the limit e/L —>■ 0, L oo {L oo first). The example which 
we will consider is a scalar field theory in 1 + 1 dimensions with polynomial 
self- interactions 



C 



m 
T 



A 
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6 J,6 



6! 



(3.30) 



In e-coordinates the Hamiltonian for this model reads after normal ordering 



dx 



4! 



(3.31) 



where the same notation as in Section |2]^ has been used. In lowest (zeroth) 
order in A4 and Ag the vacuum is the Fock vacuum, defined by afc„|0) = 0. 
This changes of course for nonvanishing couplings. For example, in second 
order perturbation theory in A4 one finds for the occupation of states in the 
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vacuum 



Po(A;,^,L) = (0|4afc|0) 



X2 1 

1111 

^ 2cj(A;) 2cj(A;2) 2uj{k3) 2uj{h) ' ' 



where = + 2em? / L, p+{q) = y—q+\Jq^ + 2em?/LjL/2e and 

= —k — k2 — k^. A similar expression is found for the (9(Ag)-term, which 
will be omitted for simplicity. In the limit L ^ oo one thus finds 




967rm^ 



(3.33) 



where 



J-cxD 277 J-cxD 2tX [Cj{z) +Uj{z2) + U){Z':i) +Cj{Zi)\ 



with Cj{z) = v^M-T and 2:4 = — z— -23- The factor L/2tt arises from going 
from discrete to continuous momentum k and we divided by the invariant 
length of the interval because the occupation in the vacuum trivially scales 
like the invariant length. 

Most importantly, the momentum scale in the occupation of the vacuum 
is set by m\/2e/L. The momentum density in the vacuum is sharply peaked 



around k = with width 0{^j2e/L) and height 0{\J L/2e), i.e. it resembles 
a (5-function as e/L 0. 

Let us now consider a state with momentum P, where P is taken inde- 
pendent of L or £. To lowest order 

\P) = aUo) (3.35) 

and thus 

pp(fc, ^, L) = {P\aiak\P) = 6k,p + O(A'). (3.36) 

Three classes of corrections contribute to pp: insertions in disconnected vac- 
uum diagrams (Fig.|3.2|a) [yielding again Eq.( p.33| )], insertions in tadpoles 
(Fig. |3.2[ b) and the rest, i.e. insertions in non-tadpole connected corrections 
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Figure 3.2: (9(A^)-corrections to the mode density in the presence of a par- 
ticle with momentum P. a.) disconnected corrections, b.) insertions into 
generahzed tadpoles (i.e. diagrams where a subgraph is connected with the 
rest of the diagram at one point only) and c.) non-tadpole connected correc- 
tions. 



The tadpole term yields 



-tadpole / 7 
Pp ['^J 



1 • ^ tadpole /I I ^ T \ 

= hm —Pp [k, \ -,L) 

L-»oo ZTT 



A4A. 



4^6 



gevrm^p 



L 



2em? 



(3.37) 



and for the non-tadpole, connected term one finds 

L 



Ppik 



lim —pfik, 



r°° dk2 
32tt J-oo 2tt 



F2 + p2 



(3.38) 



Uj{P)uj{k)uj{k2)uj{k:i) 



(k^ = P — k — ^2) plus a similar term proportional to A , which will be 
omitted in the following for simplicity. The energy denominators in Eq. (|3.38| ) , 
corresponding to the two time orderings, are given by 



Ea = P+{P) - P+{k) - p+{k2) - P+{k^) 

Eb = -p+{P)~p+{k)-p+{k2)-p+{ks). 



(3.39) 



The various contributions to the occupations in the presence of a particle 
with momentum P are shown in Fig. |3.3| for a number of values for e/L. 
The numerical values for m and P, as well as the coupling constants A4 and 
Ae, in the plots are taken to be 1. Several effects can be observed: 
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Figure 3.3: C(A^) contributions to the occupation density p{k) in the pres- 
ence of an excitation with momentum P — 1 and mass m = 1 for various 
values of the parameter e/L. Dashed hne: disconnected vacuum contri- 
bution, dotted hne: tadpole contribution, full line: non-tadpole connected 
(dispersive) contribution. 
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Figure 3.4: Same as in the previous Figure but for smaller values ole/L and 
plotted over a logarithmic momentum scale. 
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• The non-tadpole connected (dispersive) contribution scales in the limit 
L —>■ oo, e/ L ^ 0. The scaling function is the LF momentum distribu- 
tion. 

• Both, the disconnected contribution as well as the tadpole contribution, 
are restricted to a region A;^ = 0{e/L)h? near the origin, where A is 
some mass scale (A = m to lowest nontrivial order). 

• The integral over the disconnected vacuum contribution is independent 
of e/L. 

• Compared to the vacuum contribution, the tadpole term is suppressed 
by one power of \Js/ L and thus can be neglected as e/L ^ 

The gap can be most easily observed by plotting the density over a logarith- 
mic momentum scale (Fig. p.4|) . 

For very small values of e/L, the momentum distributions from the dis- 
connected diagrams [momentum scale 0{\J e/ L)] and from the dispersive 
contributions [momentum scale 0.1 — 1] no longer overlap and a gap arises. 
The disconnected contributions were already present in the vacuum (ground 
state for P = 0) and are unaltered by the presence of the excitation with 
momentum P. The only change in occupation within the small momentum 
region arises from the tadpoles but its integrated contribution vanishes as 
e/L ^ and becomes negligible in that limit. Note that y^e/L must be 
extremely small for the gap to be clearly visible. This makes nonperturba- 
tive studies of the gap forbiddingly difficult numerically because one would 
have to cover a huge number of scales (from \Js/ L to 1) while keeping the 
invariant volume large. 

3.4 Vacuum Condensates and Sum Rules 

In the previous section we explained that vacuum condensates may enter the 
effective (zero-mode free) LF Hamiltonian via induced coupling constants. 
The condensates cannot be calculated directly unless one includes dynamical 
zero-modes. However, even without zero-modes, it is possible to calculate at 
least some of the condensates indirectly using sum rule techniques. As an 
example, let us consider the two point function in a self-interacting scalar 
field theory 

G{x) = (O|0(O)0(x)|O) (3.40) 
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(x° < 0, < 0). Inserting a complete set of states one obtains 

Gix) = (O|0|O)2 + ^ r^(O|0(O)Kp)(n,p|0(a;)|O) (3.41) 

where the sum is over all particle states. The normalization of the states is 
(n, p\m, p') = 2p_(5(p_ —p'_)6nm and the energies are given by the on-shell dis- 
persion relation = M^/2p^. By boost invariance (in the continuum limit), 
the vacuum to "hadron" matrix elements are independent of the momentum 

(Ol^lrz,;.) = (3.42) 

and thus 

|2 



Gix) = (O|0|O)2 - ^ ^Ko(M„v^), (3.43) 

n 

where Kq is a modified Bessel function ||AS 70 1. In the limit — one thus 
finds 



{0\(f)''\0) - {0\(f)^\0) free = Um,2_^o [G (x) - G (x) 



free] 



= (01010)^ + E^log-^, (3.44) 

^ AtT Mfree 

where we used J2n Ifl'nP = 1 and Mfree, Gfree{x) are the invariant mass and 
the two point function for noninteracting fields. Eq. ( p.44| ) is very interesting 
because it allows to calculate (O|0^|O) in terms of (O|0|O) and quantities 
{Mn and Qn) which are calculable in a canonical LF calculation without any 
dynamical zero-modes. 

A similar trick works for the cubic condensates. Of course one has to be 
careful to separate the disconnected contributions first 

{o\<p{x)<p{y)<p{z)\o) = mo)' + mix)<Piy)<Pimh 

+ (O|0|O) [(O|0(x)0(i/)|O)c + (O|0(x)0(;^)|O)c 

+ {0\4>{y)^{z)\0)c]. (3.45) 

The connected piece is calculated similar to (O|0^|O) by inserting a complete 
set of states. In the limit (x — yY ^ 0, {y — one finds 

(010=^10) = {o\m''+mmm'\o)c 

dp 



n Jo 2p^ 

{o\m'+mmm'\o)c 

+ E^log^^, (3.46) 

n 



42 



where hn = -\/27r(O|0^|n, j9-_)c (independent of p_) and (O|0^|O)c can be 
taken from above ( |3.44|) . Note that Y.n9nhn = because the states 
J dp_ exp{ip_x~)(j){x~)\0) and J dp_ exp{ip_x~)(f)\x'')\0) are orthogonal. 
Like Qn, hn can be calculated in a LF calculation without dynamical zero- 
modes. 

The generalization of these results to higher condensates is straightfor- 
ward and by recursion one can express them in terms of (0|(/)|0) and matrix 
elements which are accessible in a LF calculation. These matrix elements [gn 
and hn) depend on the states and thus implicitly on the coupling constants 
in the effective LF Hamiltonian. Since the coupling constants in the effec- 
tive LF Hamiltonian also involve the condensates (|3.29|) , this implies that 
it may be possible to determine the coupling constants in the effective LF 
Hamiltonian self consistently. 

Similar results may be derived for Yukawa theories. In Section [4.2| we will 
relate the effective coupling constants in the LF Hamiltonian to the spectral 
densities ( [4.16| ) which are also accessible in a LF calculation. 

Extracting vacuum condensates from a canonical LF calculation via sum 
rules has for example been done in Ref. ||Zh 85|| for the ruq quark con- 
densate in QCD2(A'"c —>■ oo)- The numerical result for ijjip was confirmed 
later in Ref. ||Li 86|| in an equal time framework. A finite quark mass calcu- 
lation, based on LF wavefunctions and sum rule techniques, was first done in 
Refs. Pu 894 Pu 94b|| . Again the result agreed with the result from equal 
time quantization ||L'l' 89[|. 



43 



Chapter 4 

Perturbative Renormalization 



In practical applications of LF quantization, such as calculating parton dis- 
tributions, nonperturbative effects play a major role. Nevertheless it makes 
sense to study renormalization of LF field theories first from a perturbative 
point of view because this allows to resolve some issues which would also 
appear in a nonperturbative bound state equation. 

Most terms in the perturbation series generated by the LF Hamiltonian of 
QED or QCD are UV-divergent. This is not very surprising. After all we have 
become used to the fact that most quantum field theories contain divergences. 
However, as we will see in the following, the structure of the divergences in 
light front perturbation theory (LFPTh) is different from the divergences 
in covariant perturbation theory (CPTh). Because LF quantization is a 
noncovariant formulation of field theory, different Lorentz components of a 
divergent expression are not necessarily related to each other. In addition, in 
many examples the degree of divergence in LFPTh is worse than in CPTh. 

On the one hand this is caused by the choice of regulators. On a formal 
level, LFPTh and CPTh are equivalent ||CY 72|| . However, the "equivalence 
proof" involves steps which are ill defined in the presence of divergences 
and singularities. In practice, if one wants to demonstrate the equivalence 
between LFPTh and CPTh, it is very helpful to completely regularize the 
theory at the level of the Lagrangian — before quantizing. One possibility to 
do this is Pauli-Villars regularization, where one can introduce as many regu- 
lators as are necessary to render the theory free of divergences and light-cone 
singularities [PY 72| , |BF+ 72|| . Obviously, it is then not difficult to establish 
the equivalence between LFPTh and CPTh. However, for practical applica- 
tions, Pauli-Villars regularization is not very useful. On the one hand the 
Hamiltonian for a Pauli-Villars regularized theory is either nonhermitian or 
unbounded from below or both. On the other hand, Pauli-Villars regulators 

"'^The Pauli-Villars ghosts, must be quantized with the "wrong" commutation relations 
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are not very useful for nonabelian gauge theories, because there one would 
have to introduce massive vector fields, which will in general destroy the 
renormalizability of nonabelian gauge theories. For these reasons, one is not 
interested in employing these regulators in the context of LF quantization. 

For practical applications, it is very useful to use regulators that are 
compatible with the kinematic symmetries 0of the LF.[| In the literature one 
finds for example the Brodsky-Lepage regulator 



E^^^^<AL, (4.1) 

where the sum extends over all particles and Xi = ki_/ P^_°^ G (0; 1) are LF 
momentum fractions. Other regulators are a transverse momentum cutoff 

kl < Al (4.2) 

or dimensional regularization in the transverse direction | MP+ 91 , BL 91a|| 



Very often it is in addition necessary to introduce a cutoff for small longitu- 
dinal momenta, such as 

e(x, - 5) (4.4) 

and/or a cutoff in the number of particles (Tamm-Dancoff approximation). 
What all these regulators have in common is that they are in general not 
compatible with Lorentz transformations that are not kinematic symmetries 
of the LF (like rotations around any axis other than the 2;-axis). Thus when 
using one of these regulators, one should not be surprised if matrix elements 
do not exhibit the full Lorentz invariance — unless one compensates for this 
effect by means of a more general counterterm structure. This last point 
will be the main subject for the rest of this chapter. The Tamm-Dancoff 



approximation will be discussed in more detail in Section [5^ 



4.1 Scalar Fields 

The following observation is very helpful in analyzing the perturbative equiv- 
alence between CPTh and LFPTh: Hamiltonian (with x° or as time) per- 
turbation theory can be obtained from covariant perturbation theory after 

in order to contribute with opposite signs in loops, which is necessary to cancel the diver- 
gences. The properties of the Hamiltonian then follow from the spin statistics theorem. 

^These are all Poincare transformations, which leave the x"*" = initial surface invariant, 
such as translations, rotations around the z — axis or longitudinal boosts. 

■^This excludes, e.g. Euclidean lattices. 
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Figure 4.1: 1-loop self-energy diagram in 0f+i. 

integrating the energies (i.e. the momentum variable which is canonically 
conjugate to the "time") first. Thus from the mathematical point of view, 
the question about equivalence between LFPTh and CPTh has been reduced 
to the question whether the order of integration plays a role in a Feynman 
integral. 

As an example, let us consider the 1-loop self-energy E in {/(^-theory in 
1+1 dimensions (Figure |0| ) 

(4.5) 
(4.6) 



ig f dk_dk+ 1 



2 J (27r)2 k"^ — m'^ + ie {p — k)'^ — m? + ie 

r dk_ Q{k_) e(p_ - k_) 1 

~2 J 2tt 2A;_ 2(p_ - A;_) p+ - ^ ^ 



2k- 2{p--k-) 

First, without going into the details, it is easy to convince oneself that 
Eq. (|4.6|) is exactly what one obtains in LF-Hamiltonian perturbation theory: 



P+ ~ li 2{p^-k ) energy denominator and the 9-functions ensure 

that all momenta are positive. The other factors arise from a vertex factor 
proportional to {k- (p_ — k_)) ' at each vertex. It is also easy to see that 
Eq. (|4.6|) agrees with the covariant calculation with symmetric integration. 
After substituting /c_ = xp_ in Eq. ([4.6|) one finds 



S = ^r^ ^ . (4.7) 

2 Jo An p^x{l - x) - ^ ' 

In the covariant calculation one first combines the two denominators in 
Eq.( |4.5|) with a Feynman parameter integral and then one integrates sym- 
metrically over d'^k. This reproduces Eq.( [4.7D where the x-integration corre- 
sponds to the parameter integral. 

Our next example will be one where the order of integration does matter, 
namely the so called simple tadpole diagram in cf)^ (for simplicity again in 
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Figure 4.2: Typical generalized tadpole diagrams for 0^ 



1+1 dimensions) 

^ ^9 f d'k f 1 1 \ .4o^ 

2 7 (27r)2 U2 - m2 + p-Aa + .^y' ^ ' 

We have already performed a subtraction because the unregularized in- 
tegral diverges logarithmically. Symmetric integration over d'^k yields 
S = ((yf/87r) log A^/m^. In LFPTh (unsymmetric integration; fc+-integral 
first) one obtains zero: for k_ ^ one can always close a contour inte- 
gral in the complex /c+ plane such that no poles are enclosed. The surface 
term vanishes because of the subtraction term. The point A;_ = is usually 
omitted in LP quantization without zero-modes. The mathematical reason 
for the difference between the LFPTh result and the CPTh result is a term 
oc 6{k^), which is omitted if one (as is usually, either explicitly or implicitly 
done) has a small k_ cutoff, like B(A;_ — e) at each line — even in the limit 
e^O jCM 6911 . 

This result is very typical for pathologies of LFPTh with scalar fields. 
Compared to CPTh, one omits certain diagrams which are nonzero in CPTh, 
i.e. LFPTh yields a priori wrong results! Fortunately (later we will see 
that there is a good reason for this) the 'mistake' does not depend on the 
external momenta. Thus one can make up for the mistake by means of a 
local counterterm in the Lagrangian. Other diagrams which suffer from the 
same problem are the generalized tadpole diagrams, i.e. diagrams where part 
of the diagram is connected with the rest of the diagram only at one single 
point (examples are shown in Figs. [4.2| and ^.3| ). As discussed in Ref. |[Bu 93 



they are all zero in LFPTh. However, because the generalized tadpoles in 
these diagrams are connected to the rest of the diagram only at one point, the 
covariant calculation yields a momentum independent result for the tadpole 
part (just a number), which can thus always be replaced by a local insertion 
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Figure 4.3: Typical tadpole diagrams arising for scalar fields with more gen- 
eral polynomial self-interactions 



into the diagram. In practice this means that the fact that all generalized 
tadpoles are (wrongfully) zero on the LF, can be easily compensated by 
appropriate redefinitions of bare coupling constants! Furthermore, tadpole 
diagrams are the only diagrams which are treated incorrectly in naive LFPTh. 
A very interesting result is the relation between the tadpole counterterms 



and vacuum condensates pu 93|| . For example, each tadpole correction to 
the propagator in Fig. [4.2| a can be written as a mass insertion times the free 
field vacuum expectation value (VEV) of (O|0^|O). The generalized tadpole 
in Fig. [4.2| b corresponds to a mass insertion times a higher order correction 
(O|0^|O). The higher order tadpoles in Fig.|i7^ correspond to mass (a) and 
vertex (b) insertions times a term that contributes to (Ol^'^lO). Suppose the 
interaction term in the original Lagrangian is 

A„* = -EV- (4-9) 

n "'■ 

Then all the 'missing tadpoles' are automatically taken into account if one 
uses 

In other words, Cint,eff with naive LFPTh yields the same results as Cint 
( [4.9|) ([4.101) with CPTh to all orders in perturbation theory, if the VEVs in 
Eq. (|4.10|) are also given as a perturbative expansion (calculated in CPTh) 

First of all this result is very useful in practice, because, given the orig- 
inal interaction, it allows one immediately to write down an ansatz for the 
effective LF interaction — even if the VEVs cannot, in general, be calcu- 
lated from the LF Hamiltonian. Secondly, although derived perturbatively, 
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Eq. (|4.10|) formally agrees with Eg. ( p.29| ), which was derived nonperturba- 
tively using e-coordinates. Of course, while Eq.( 4.10| ) was derived only for 
cases where the VEVs can be calculated perturbatively, Eq.(|3.29|) is valid 



in general. However, the formal agreement between the two results gives us 
confidence to approach other zero-mode problems using perturbation theory 
as well. 

As it stands. Eg. ( 4. 10 ) is valid only for superrenormalizable theories be- 
cause we have only addressed longitudinal divergences in the above discus- 
sion. For renormalizable theories one must first cut off the transverse diver- 
gences, e.g. by using a transverse lattice (see Section |5.3|) or dimensional 



regularization in the transverse direction ||MP+ Qlj [BL 91a|| . However, with 



such a transverse cutoff in place, Eq.( |4.10| ) is valid for field theories in more 
than one spatial dimension as well. 



4.2 Fermions 

For the applications of LF quantization to DIS, we are of course not inter- 
ested in self-interacting scalar fields but rather in theories with fermions and 
gauge fields. As a first step towards this direction, let us consider fermions 
interacting with pseudoscalar mesons via a Yukawa coupling (see also Section 

Cint = 9p^il5^X (4.11) 

within the framework of LFPTh. First one may be tempted to expect that 
abovementioned perturbative zero-mode problem does not occur here, be- 
cause a priori there are no tadpoles in Yukawa theory with 75-coupling. 
However, after eliminating the non-dynamical component of the fermion field 
from the theory, the canonical LF-Hamiltonian (|2.28|) does contain 



terms which are fourth order in the fields — giving rise so-called "seagull" - 
diagrams (Fig. [4.4|) . It is thus not very surprising that the perturbative zero- 
mode problem arises in diagrams which have the topology of a seagull with 
one vacuum contraction. (Figure because this is the topology one ob- 



tains if one replaces eithe r ipj^'^d_ ^tp(^^) or 0^ by their VEVs. In practice, 
this works out as follows ||BL 91a| , [BL 91b| , [LB 93|| : consider, for example. 



the dressed one loop self-energy diagram for a fermion 

(4.12) 



^Here we assume self-consistently that all sub-loop countcrtcrms have been added to 
the LF result, such that the full fermion propagator is covariant. 
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Figure 4.4: Seagull diagram in x'^ ordered perturbation theory, representing 
four-point interactions induced by eliminating V'(-)- The dashed lines are 
bosons and the full lines represent fermions. The "slashed" fermion line 
corresponds to instantaneous (with respect to LF-time) fermion exchange. 



Figure 4.5: Typical self-energy diagrams in Yukawa theory, which have 
the same topology as a seagull with one contraction. The blob stands for 
arbitrary self energy insertions. 
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where the spectral functions (p, pi, P2) parameterize the (unspecified) self- 
energy insertions. They satisfy (follows from the canonical commutation 
relations) /q°° (ip^p(p^) = dm^ piirn?) = 1. As far as the k+ integral is 
concerned, the most singular term in S is the one proportional to 7^fc+. We 
thus consider ^ 



tr(S7 



9l 



d'^k 



(2vr) 



D 



(4.13) 



dfi'^ 



p{p^ 



{p — ky — p^ + ie Jo 



dm^ 



k+pi{m^) 

f^2 _ ^2 _|_ 



To identify the troublemaker we eliminate k^ in the numerator using the 
algebraic identity 



k. 



{k"^ — w? + ie) {{p — ky — p2 + ie) 



2(p_ — kJ)p^ — 'P± + 2p_L ■ A;_L + — p^ 



2p_ 



2p_ 



{k"^ — rn? + ie) {{p — k)"^ — p^ + ie) 
1 1 



((p — A;)2 — p^ + ie) {k"^ — w? -\- ie) 



.(4.14) 



The important point here is that the last two terms in Eq. ( [4.14| ) give 
(5-functions in p_ — /c_ and A;_ respectively after the /c+ integration. These 
5-functions are missed in the naive LF Hamiltonian without zero-modes (this 
is very similar to the tadpoles in self- interacting scalar fields). One finds (we 
subtract here the one-loop result because this allows to drop the surface term 
in the complex fc_ plane; pf^'^^{p'^) = 5{p^ — pi), p('"'^^(m^) = (5(m^ — ml)) 



■loop 



covariant 



9l 



d^k 



loop 



canonical LF 



J {27T 







0-p/-^(p2) 



9% 



2p_ J {2tt 



d^k 



k'^ — p^ + ie 

2\ 



pfrec(^^2^ 



k^ 



+ ie 



(4.15) 



Since the other component of S have no problems from zero-modes this 
immediately implies 



^1 ^^^^v) covariant ^1 



loop ) canonical LF 



more detailed study shows that the other components are free of trouble [BL 91a 



BL 91b 
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g%,- f d-k r,^.PiK)-pf^K).(4.i6) 



2p„ y (27r)4yo A;2-m2 + ie 

This result is very interesting for the following reasons: 

• canonical LF quantization disagrees with covariant perturbation theory 

• the mistake of canonical LF quantization can be compensated by a 
counterterm to the mass term in the kinetic energy (but not the mass 
term appearing in the vertex) 

• if one adds the wrong counterterm, rotational invariance and parity 
invariance for physical observables are broken. This can be used as a 
renormalization condition to "fine-tune" the coefficient of the countert- 
erm. 

• the counterterm is related in a simple way to the spectral function of 
fermions and bosons which are numerically calculable in a canonical 
LF-calculation! 

• The boson contribution in Eq.( [4.16D can even be expressed in terms 
of a local VEV: = (0| : 02 - |o)^2 7+/2j9_. Unfortunately this 
is not possible for the term containing the fermionic spectral density, 
which would read 511^^^" = (0| : ip^^p : |0)^|,7+/2p_. 

Note that, in order to obtain the full counterterm necessary to establish 
agreement between a covariant calculation and a canonical LF calculation, 
one still has to add the one loop counterterm — but this should be obvi- 
ous and can be easily done. Similar statements hold for fermion loops in 
the boson self energy. The only difference to the above example is that the 
difference between a covariant calculation and a canonical LF calculation re- 
sults in a difference in the bare boson mass; i.e. no space time symmetries 
can be used to fine-tune the counterterm. However, the difference can still 
be related to the spectral density of the fermions. Besides the "contracted 
seagulls" , only disconnected vacuum diagrams — which are irrelevant for the 
dynamics of physical states — suffer from the zero-mode problem. It is thus 
also sufficient to tune the vertex mass and the kinetic mass independently 
and those masses and the boson mass independently from the corresponding 
coefficients in the covariant Lagrangian in order to recover equivalence be- 
tween covariant calculations and canonical LF calculations. Note however. 
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that (like in the self-interacting scalar theory) all this holds only after render- 
ing the transverse momentum integrals finite (e.g. by means of dimensional 
regularization ||MP+ 9l| , |BL 91a|| or a transverse lattice (Section |5.3| ). 

It should also be noted that perturbative zero-modes also play a role in 
higher-twist parton distributions. There they can lead to violations of naive 
sum rules as discussed in Refs. Pu 92b|, Pu 92c|, Pu 95||. 



4.3 Gauge Theories 



In gauge theories the situation is much less clear than in scalar field theories or 
Yukawa theories, because of notorious infrared singularities in the LF gauge 
A'^ = 0. Certain attempts have been made to perturbatively renormalize 
LF-QED ||iVlP+ 9l| , [BL 914 [EFTi^ and QCD PZT)^ ]. 



In the context of calculations of the electron's anomalous magnetic mo- 
ment in QED it has been shown (up to three loops in Feynman gauge and 
up to two loops in LF gauge) that all A;_ ^ singularities in LFPTh cancel 
— provided one adds up all diagrams that contribute to a given order in the 



coupling constant [|LB 93|| . The regulators used were Pauli-Villars regulators 
or dimensional regularization in the transverse direction. Furthermore only 
two extra Q counterterms are necessary to render the theory UV- finite: a 
kinetic mass counterterm for the electron (similar to the one discussed in 
Section |4.2| ) and a mass term for the transverse photon field. The numerical 
result for {g — 2) thus obtained agrees with the known result from covariant 
calculations. 

Perturbative LF calculations of vertex functions, which employ a Tamm- 
Dancoff truncation were done in Ref. ||MP+ 9i| for QED and in Ref. ||HZ 93 



for QCD. Due to incomplete cancelations of A;_ ^ singularities in the 
Tamm-Dancoff approximation, infrared singular counterterm functions were 
already in lowest nontrivial order necessary to render the results finite. 

The first calculation, relevant for asymptotic freedom, was performed 
in Ref. ||'Lh 79|| (four gluon vertex) and Ref. ||CI''+ 80|| (quark-gluon ver- 
tex). Further discussions on renormalization on QCD in LF gauge (but not 
LF quantization) can be found in Ref. | BD+ 85|| - For demonstrations of 
asymtotic freedom, employing both LF gauge and LF quantization, see Ref. 
Pe 94b|| (and references therein). 



In both types of calculation (LFPTh and LFTD) even perturbatively the 
structure of the renormalized LF-Hamiltonian is not known to higher orders. 
Perhaps the cleanest way to address the problem would be to start from the 
axial gauge in e-coordinates in a finite box ||LN+ 94a] , [LN+ 94b| , [LN+ 94c 



^That is beyond those counterterms which are required in a covariant calculation. 
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and to approach the LF by carefully taking the limit e — (as in Section 
3.3). While it is conceivable that this is feasible in QED, the axial gauge 
Hamiltonian for QCD^^i in a finite box ||LN+ 94^] is perhaps too complicated 
to allow one to study this limit with appropriate care. 



4.4 Summary 

In the renormalization of LF field theory, one can distinguish three kinds 
of counterterms. First the usual renormalizations, which can be handled by 
making the bare coupling constants in the Lagrangian cutoff dependent. In 
the following these will be referred to as canonical counterterms. Secondly, 
counterterms that have to be added when one is employing a Tamm-Dancoff 
cutoff. These will be discussed in Section |5.4|. Typically, one needs an infinite 



number of counterterms! Third, effects caused by an improper treatment of 
zero-modes in the canonical approach. In those cases were these effects are 
now understood the renormalization of zero-mode effects can be accomplished 
by adding a finite number of counterterms that have the structure of tadpole 
and seagull diagrams with some lines "ending in the vacuum" . In general, the 
zero-mode counterterms are already included in the list of "Tamm-Dancoff 
approximation counterterms". This means zero-mode effects become irrele- 
vant when one uses a Tamm-Dancoff approximation. However, in the absence 
of a Tamm-Dancoff approximation, i.e. in calculations without or with negli- 
gible restrictions on the Fock space (see Section |5.3| ) or in perturbation theory 
if one adds all diagrams to a given order in the coupling, abovementioned 
tadpole or seagull counterterms are quite relevant because they are the only 
counterterms needed besides the canonical counterterms. 
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Chapter 5 

Nonperturbative Calculations 



5.1 Discrete Light-Cone Quantization 

The most straightforward method for solving bound state problems in the 
context of LF quantization is discrete light-cone quantization Q (DLCQ) 
[pB 85|| . For extensive reviews and more references see Refs. [|BP 91|, [Pa 93|, 



The basic idea in DLCQ is as follows (for simplicity we illustrate the 
method using the example of One puts the system into an x~-box of 

length L with periodic or antiperiodic boundary conditions 

0(x" + = ±0(x",x+). (5.1) 

In the following, antiperiodic boundary conditions will be used, which implies 
for the mode expansion 



[x 



(5.2) 



= ^ f - 1 ) . (5.3) 



where 

L \ 2, 

The main reason for choosing antiperiodic boundary conditions is that 
one does not have to worry about the mode with p_ = 0. Another reason 

^Like the canonical quantization discussed in Chapter H, the quantization surface in 
DLCQ is the plane x'^ = 0, i.e. a front or plane — and not a cone. Thus discrete light 
front quantization (DLFQ) would be a more appropriate terminology. However, because 
of historical reasons, the method has been named DLCQ in the literature. 

^In the presence of interactions which contain odd powers of (j) one has no choice and one 
must use periodic boundary conditions — otherwise momentum conservation is violated 
at the boundary! 
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is that many numerical problems converge faster when antiperiodic bound- 
ary conditions are used (compared to periodic boundary conditions with the 
p_ = mode left out). This can be understood in perturbation theory be- 
cause there are often non-negligible contributions to Feynman integrals from 
the region near p„ = 0. Let / be some typical function that appears as the 
argument of some Feynman integral. Then eJ2'^=-oo f ((^ ~ |)^) usually 
a better approximation to /(x) than e Y^~t-oo f (^^) + ^ J2'^=i f (^^) 
because in the latter expression the point n = is missing compared to the 
trapezoidal quadrature formula. 

In order for 0(x) to satisfy the canonical commutation relations (see Chapter 
0)' 

[d.<P{x), mU^y^ = -'-6{x- - y-) (5.4) 
we impose the usual commutation relations for the coefficients a^, 

5kq. (5.5) 



The above expansion is then inserted into the momentum operator 

P_ = /" dx~ : d^chd^d) : 





— - - j (5.6) 



and the Hamiltonian 







where 



is the kinetic term and 



^(T + V), (5.7) 



^2 oo t„ 
^ fc=l 2 



8^4! kuk,MM=i ^Jh-l ^h-\ 

(5.9) 

is the interaction term. 5pJP^ is a momentum conserving Kronecker 5. Since 
the length of the box completely factorizes, it is useful to work with the 
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rescaled operators 



K 
H 



(5.10) 
(5.11) 



Since the momenta of all excitations are discrete and positive, the Fock space 
is finite dimensional for all K. Thus, at least in principle, one can now pro- 
ceed as follows: for fixed K {K and H commute) one diagonalizes H (which 
is a finite matrix for finite K). From the eigenvalues Ei one computes the 
invariant masses = 2KEi and from the eigenstates one can compute 
other physical observables (like parton distributions). In general, physical 
observables thus computed will of course depend on the "resolution" K. The 
continuum limit is obtained by extrapolating to — > cxd. The diagonaliza- 
tion is generally done using brute force matrix diagonalization or, if one is 
only interested in the lowest states, using the Lanczos algorithm |pi 91 . 



At this point one encounters a problem that is inherent to Hamiltonian 
systems: the dimension of multi-particle states in the Fock space expansion 
grows exponentially with the number of particles. The number of particles, 
as well as the number of states for a single particle are both limited by 
the longitudinal momentum K, i.e. the dimension of the Fock space basis 
shows factorial growth with K. Fortunately, in 1 + 1 dimensional exam- 
ples, the factorial growth sets in only rather slowly and numerical conver- 
gence for typical observables can be obtained before the size of the matrices 
becomes a problem. DLCQ was enormously successful in many 1 + 1 di- 
mensional field theories |[EP+ 87[ [HV 87| , |EP 891 |Nfc88| , |Bu 89b| , PP+ 90 



|Bu 89a] , [Bu 93| , |Bu 89c| , [BB 91|| . In all cases, where results from other ap- 
proaches to field theories were available agreement could be shown within 
numerical uncertainties (QEDi+i: ||(Jo 76|| vs. ||EP+ 87| , |EP 89|| Q ,QCDi+i 
|Ha 82|| 



vs. 



|HP+ 901 , |Bu 89a|| , sine-Gordon model: ||DH+ 75 



vs: 



Beyond reproducing known results, DLCQ has been used to calculate new 
and interesting results in QCD]^_,_]^: the most notable results are the exis- 
tence of a nucleon-nucleon bound state and the analysis of the nuclear quark 
distribution in comparison with the nucleon quark distribution. Not only 
exhibits the 1 + 1 dimensional "deuteron" an EMC-effect, but it can also 
be understood analytically due to the simplified dynamics in 1 + 1 dimen- 



sions 



Bu 89b|| . Typical Euclidean lattice calculations are too "noisy" to even 



■^However, there is still a 1% difference in the fundam ental m eson mass for the term 
linear in nie in QEDj^^j^ as calculated from bosonization |Co_7£] and in LF-quantization 



I Be 77 1 ■ It is not clear whether this deviation is due to the finite Fock space truncation or 
whether this is a real problem Ma 92, 5t 92 . 
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demonstrate binding of hadrons. Another remarkable result from DLCQ cal- 
culations in QCD^^^ dimensions is "Anti-Pauli-Blocking" [ Pu 89c| , PB 91|] : 
contrary to the naive expectation, sea quarks in nucleons in QCD^^^ tend to 
have the same flavor as the majority flavor among the valence quarks (i.e. 
more u than d in a nucleon ipvaience = uud). 

In 2 + 1 or 3 + 1 dimensions the situation changes drastically, because 
there the exponential growth is much more rapid. The basic reason is that 
there are now transverse degrees of freedom besides the longitudinal degrees 
of freedom. Suppose that each particle can occupy N states for each spatial 
dimension. Then the Fock space basis size grows like N^^p"-^* with the num- 
ber of particles in 3 + 1 dimensions, while the corresponding growth would 
be only N^p"'"^ in 1 + 1 dimensions. For a concrete example (0^ with an- 
tiperiodic boundary conditions in the longitudinal direction) this works out 
as follows. For a longitudinal momentum K = ^ longitudinal momen- 
tum states accessible) the Fock space basis size is 27. If one has just two 
transverse degrees of freedom (e.g. two points in the transverse direction) 
the basis size grows to 426. For 8 x 8 = 64 degrees of freedom in the trans- 
verse direction, that number grows to 6 ■ 10^^. These astronomical numbers 
clearly demonstrate that any direct matrix diagonalization approach or even 
a Lanczos type algorithm is doomed to fail because one is not even able to 
store the wavefunction in a computer ||Bu 94a|| . 

The most simple (and perhaps most drastic) way out of this dilemma is 
to impose additional cutoffs, like restricting the number of particles. Typ- 
ically, this means restricting the Fock space to 3 (perhaps 4) particles or 
less [ [KP 92| , |Wo 93|| . In QED, since the coupling is small, this is a good 
approximation. However, to the same order in a within which the 3 particle 
truncation is a good approximation one can calculate the parton distribu- 
tions analytically ||Bu 92a|| . That is, even in QED there is not much point in 
doing numerical DLCQ calculations with Fock space truncations to the low- 
est nontrivial order! In QCD, where one faces an intrinsically strong coupling 
problem, restricting the Fock space to the lowest nontrivial component seems 
entirely useless. For example, even if one allows up to 4 particles (which is 
about the maximum that can be handled numerically using the Lanczos al- 
gorithm), this means one allows at most one gluon in addition to the three 
valence quarks in a proton. That is there is no chance one can "see" any 
effects from nonlinear gluon-gluon couplings. 

It should be emphasized, that this problem is not specific for LF field the- 
ories, but occurs in many Hamiltonian approaches to field theory — and in 
many cases could be solved. Thus there are many numerical methods avail- 



A caveat to this pessimistic point of view will be discussed in Section 5.4 
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able which can potentially be useful in overcoming the difficulties associated 
with exponential basis size growth. 



5.2 Functional Integration on a 
Longitudinal Lattice 

Functional integrals on Euclidean lattices have been very successful in solving 
ground state properties of QCD (e.g. vacuum properties, hadron masses and 
ground state matrix elements). However, since two points on a Euclidean 
lattice are always separated by a space-like distance, it is only very indirectly 
possible to extract information about light-cone correlation functions from 
these calculations. Of course, this is because in conventional Euclidean field 
theory exp(— /?P^) is used to project on the ground state wave function 
of P° at equal time. Thus as a caveat one might be tempted to consider 
a similar formalism for LP Hamiltonians. Suppose one discretizes the x~ 
direction, | and uses a functional integral to project on the ground state of P+ 
Mu 8^ , |HL 93|| . This results in an immediate problem because the LF-energy 



decreases with increasing momentum (P+ = M^/2P_ in the continuum, on 
a lattice there is a minimum for P_ = 0{l/a)). Due to Bragg reflections, 
momentum is not conserved and the particles tend to accumulate near the 
minimum. However, since the momentum near that minimum is of the order 
of the inverse lattice spacing, the particles always "see" the lattice and no 
meaningful continuum limit is obtained. It is conceivable that this problem 
can be cured by adding a Lagrange multiplier proportional to the total LF 
momentum to the lattice action (in the continuum limit this amounts to 
minimizing P+ = P+ + AP_ instead of P+). However, this idea will not 
be investigated here any further. Another difficulty of the longitudinal LF- 
lattice is species doubling for bosons |[Mu 88||! 



5.3 Hamiltonian Monte Carlo on a Trans- 
verse Lattice 

While Mont Carlo calculations for longitudinal LF-lattices seem to be 
plagued with difficulties, this is not the case for the transverse lattice 
[ Pu 94a|| . On a transverse lattice one keeps the longitudinal directions (a;+ 



and X ) continuous, while discretizing the transverse coordinate (Fig. |5T 
BP 76|| For simplicity, let us consider self- interacting scalar fields in 2-1-1 



^The transverse coordinates are irrelevant in this argument. 



59 



J- space 
(discrete) 



long, space 
(continuous) 



time 
(continuous) 



Figure 5.1: Space time view of a transverse lattice 



dimensions on such a transverse lattice characterized by the action 



A 



cont. 



(fx 



fi=0 



2 ^ 



int 



(0) 



(5.12) 



Upon discretizing the transverse direction (spacing a) one thus obtains 



A' 



latt. 



d^x 



1 

E 

ti=0 



d^Cpnd^Cpn (0n+l " 0n)' 



2a2 



m 
T 



Hi 



int \ 



(5.13) 

Up to a factor of a (which can be absorbed into a redefinition of the field 0„), 
Eq.( |5.13| ) looks like the action for a multi-flavor theory in 1 + 1 dimensions 
(where the site index n corresponds to the "flavor" index). In the next 
step one constructs the DLCQ Hamiltonian for this "multi-flavor" 1-1-1 
dimensional theory. The important point here is that the action is local in 
the transverse direction, i.e., there are only nearest neighbor interactions. 
Since the DLCQ-Hamiltonian is thus also local. 



(5.14) 



one can apply many Monte Carlo techniques which have been developed for 
other Hamiltonian systems (see e.g. Ref. | N0 87 |). One technique which 
turns out to be particularly useful for LF-Hamiltonians on a transverse lat- 
tice is the ensemble projector Monte Carlo technique [PP 85|| based on the 
so called checkerboard decomposition of the Hamiltonian [|HS-|- 82||. Using 
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locality of the Hamiltonian one can write 



HDLCQ = Ha + H, (5.15) 

where 

Ha = X! 51 = -f^l,2 + -f^3,4 + ••• 

n=l,2,3,... ^ n=l,3,5,... 

= E ^+ E K,„+l = i^2,3 + i^4,5 + - (5.16) 
n=l,2,3,... ^ n=2,4,6,... 

In other words, the DLCQ Hamiltonian can be written as a sum of two 
terms, each of which can be written as a direct sum of two-site-Hamihonians. 
The point to all this is that while the dimension of the space on which 
Hdlcq acts is astronomical, the two-site-Hamiltonians act only on a very 
small Hilbert space (for our above example with K = ^ and say 16 transverse 
sites: dim{Hj:,icQ) = 7.8 ■ 10* but dim{Hn^n+i) = 426). The method is 
called checkerboard algorithm because one approximates the time evolution 
operator of the system by alternating infinitesimal time evolution operators 
generated by Ha and Hb respectively 

^-eHoLCQ = Q-^^^e-'^'e'^"^ + C(e^). (5.17) 

If one axis of the checkerboard is the discretized space direction and the other 
the time, Eq. (|5.16|) can be interpreted as if interactions between sites occur 



only across the black squares |[HS+ 82| . 



Before we explain how the infinitesimal time evolution operators are mul- 
tiplied together, let us pause here for a moment and understand the ad- 
vantage of the transverse lattice with DLCQ over the longitudinal lattice 
discussed in the previous Section. The main cause of the problem in the pre- 
vious Section was lack of longitudinal momentum conservation. In DLCQ 
P_, the longitudinal momentum, is manifestly conserved. Furthermore, P_ 
is just the sum of momenta at each site 

P_ = E^-, (5-18) 

n 

and the checkerboard algorithm is compatible with with longitudinal momen- 
tum conservation {[Ha, P-] = [Hb, P_] = 0). I.e. with DLCQ on a transverse 
lattice, longitudinal momentum is conserved at each step of the calculation. 
Hence, one can minimize P+ while keeping P_ manifestly fixed and there are 
no "runaway solutions". 

In the actual calculations one uses Monte Carlo techniques to calculate 
f^-£H\ \iIj-(^K)) by alternate application of e~^^°-, e~^^^, e~^^",...,e~§-'^" to 
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\%lji[K)). Here \ipi{K)) is an initial guess for the ground state wave function 
with longitudinal momentum K. For N oo one thus obtains an approxi- 
mation (because e is finite, the result is not exact) to the ground state hadron 
with the same good quantum numbers as \il)i{K)). One very useful tech- 
nique is the ensemble projector Monte Carlo method ||1NU 87| , pP 85[| , which 
works as follows for these systems: 

1) let |n) be a complete set of states 

(here product basis of Fock state bases at each site) 

2) make a good guess for \ipi{K)) 

(here a valence state with Pj_ = : \ilJi{K)) = J2nj_ (^^K,n^\^)) 

3) start from ensemble of states (from set |n)) 

4) for each select a new state li^J^) with probability 

and calculate the score 

1^(^«,^)) 

note: W{i^J-\i^^^) factorizes into two-site probabilities 
local (one pair of sites at a time) "updating" possible 

5) replicate states with multiplicity: 

-,(k,k-l) 



mt 



— + random number G (0; Ij 



s 

{S: av. score) 

Thus paths with large scores contribute with multiple weight, while 
paths with small scores get eliminated. 

6) repeat while alternating Ha and Hh 
(| in exponent only in 1** and last step!) 



^Of course, only those quantum numbers which are associated with exact symmetries 
of DLCQ on a transverse lattice (like C-parity or baryon number in QCD) are relevant 
here. 
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7) observables from ensemble average, e.g., energy of ground state hadron: 



Eq = lim 



Y.Mf{K)\H\^^^))sign 


'c{N,N-l) . 


. c(i,o)" 




■g(N,N-l) _ 







other observables (e.g. an observable diagonal in the basis) 

{MK)\d\MK)) = 



lim 

M,N^oo 





' q{N+M,N+M-1) . . ^(1,0)' 




6 




Eu{MK)\i''''''^)st9n 


' QiN+M,N+M-l) r.(l,0) 





In this Monte Carlo procedure, one only has to store the ensemble of states 
at one "timeslice" plus the result of the measurement of the observable after 
slices. Thus, at least in principle, one can handle very large lattices. The 



main advantages of the transverse lattice are as follows pu 94a 



longitudinal momentum is manifestly conserved no runaway solu- 
tions 

parton distributions are diagonal in the DLCQ-basis 

LF-vacuum is trivial — no statistical fluctuations from updating the 
vacuum far away from physical states on huge lattices. 

species doubling for fermions occurs only for the latticized transverse 
dimensions can be easily compensated by staggering [|Gr 93||. 



• excited states are suppressed by the square of their masses: 
exp(-A^£P+) = En 1^) exp{-NeM^/2P_){n\ instead of 
exp(— A^eP^) = Y.n l'^) exp(— A^£:M„)(n| which one encounters in a con- 
ventional Hamiltonian formulation. 

It is interesting to see how confinement emerges on the transverse lattice 
in the limit of large lattice spacing: In this limit, the coupling between the 
sheets is weak and the energy scale associated with link field excitations is 
high. Therefore, when one separates two test charges in the longitudinal 
direction, the transverse lattice behaves similar to QCDi+i and linear con- 
finement results trivially. For transverse separations between the charges, 
a different mechanism is at work. Gauge invariance demands that the two 
charges are connected by a string of link fields. In the limit of large spacing 
the link fields fluctuate only little and the energy of such a configuration can 
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be estimated by counting the number of link fields needed to connect the 
charges, which again yields linear confinement. 

Some of the disadvantages of the transverse lattice are: Since — > ix~^ is 
not a Wick rotation (it is just a mathematical trick to project on the ground 
state of -P+), the metric is not Euclidean and thus propagators oscillate. 
Hence, negative scores occur already for bosons which leads to an increase 
in the statistical fluctuations. However, these negative scores turn out to 
have only a small statistical weight and the resulting "sign-problem" is not 
serious. Very often in LF calculations, large cancellations occur between 
different terms in the Hamiltonian. For example, the instantaneous photon 
exchange has a 1/g^ singularity which is canceled by vertex factors in photon 
exchange. In general, it is difficult to obtain such cancellations from a Monte 
Carlo calculation. Another difficulty is that gauge invariance on a lattice 
can only be maintained if one introduces link fields. On a transverse lattice 
this amounts to introducing 1 + 1-dimensional gauged nonlinear sigma model 
fields on each link ||BP 76|| . Constructing a Fock space basis out of these 



nonlinear degrees of freedom and calculating appropriate matrix elements is 
a nontrivial task [pr 92b| , pr 94b |. 

The sign problem associated with fermions is a notorious difficulty for 
Monte Carlo algorithms: due to the minus sign in exchange terms, the in- 
finitesimal time evolution operator tends to contain many negative matrix 
elements. This very general problem is also expected to afflict Mont Carlo 
calculations on transverse lattices. However, since the LF vacuum is trivial, 
there are no sign fluctuations from Z-graphs and vacuum diagrams. Thus one 
expects that the sign problem on the LF is less severe than usual. Whether 
this improvement is sufficient to render fermions tractable on transverse lat- 
tices has not yet been investigated. 

Obviously, the transverse lattice lacks manifest rotational invariance 
which must be restored in the process of renormalization. Recently, a tech- 
nique has been described that allows easy computation of the potential be- 
tween infinitely heavy quarks in a LF framework |Pu 94c|| . Demanding rota- 
tional invariance for this observable may prove to be a powerful tool in such 
a procedure. 



5.4 Light-Front Tamm-Dancoff 



As we have discussed already in Section |5.1| , the dimensionality of the Fock 
space grows dramatically as one includes higher Fock components. Clearly, 
since as is fairly large at a low momentum scale, a numerical solution of 
bound state problems in QCD (which includes all scales) necessarily involves 
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many Fock components. In this chapter we will discuss the light-front Tamm- 
Dancoff (LFTD) approach to LP problems (for a comprehensive review see 
Ref. ||WW+ 91 ). The basic idea is very simple ||PH+ 90| , |WK 94|| : Hadrons 
are complicated objects only if one tries to build them in terms of bare quarks 
and gluons whose masses and couplings are renormalized at a scale of 1 GeV 
or higher. In terms of collective excitations (constituent quarks) ground state 
hadrons are rather simple. One of the problems with the constituent quark 
model is that a priori the interactions among the quarks are ad hoc. 

The goal of LFTD field theory is to systematically eliminate higher Fock 
components and high energy degrees of freedom[]. As one goes to lower and 
lower scales the interaction between the (dressed) constituents thus becomes 
more and more complicated. If the whole program is successful, constituent 
quarks will emerge as the quasiparticles of QCD at intermediate energies. 
A major virtue of using LF quantization in this approach is that it stays 
close both to physical intuition (which may prove very helpful when it comes 
to developing variational methods to analyse the Hamiltonian) as well as to 
experimental observables at large momentum transfer (useful for phenomeno- 
logical applications). 

A systematic Fock space expansion, based on a Hamiltonian formulation, 
for field theory was originally developed by Tamm [P?a 45|| and independently 
by Dancoff [Pa 5CI|| . It turns out that such an approach is doomed to fail if the 
perturbative ground state (the starting point of the expansion) is too far from 
the actual ground state. In such a situation one needs very (or infinitely) 
complicated Fock states just to build the ground state. LF quantization 
is advantageous at this point, because the vacuum of a LF Hamiltonian is 
trivial. 

In fact, because of the vacuum, LF quantization is probably the only 
framework, where such a programm can possibly work. 

In practice, even within LF quantization, it is of course not possible to 
integrate out high energy states and higher Fock states exactly. Instead one 
writes down a catalog of all interaction terms that are allowed by power 
counting [|Wi 91| , pP 92|| : on the LF there are two length scales x~ and x±. 
The engineering dimensions of the various operators and terms that enter 
the LF-Hamiltonian in 3 + 1 dimensions can be easily derived from free field 
theory ||GP 92|| (0 stands for a scalar field or for A± — the transverse gauge 
field components which have the same engineering dimension as scalar fields; 



^This procedure is e xplicite ly demonstrated for the simple example of in the two 
particle sector in Ref. [ Am 94 1. 
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ip^^^ is the dynamical component of a fermion field) 



■ 1 ■ 


, 5± = 


■ 1 ■ 




■ 1 ■ 






, m = 




.x~ . 











{+) 
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, = 







The Hamiltonian and the Hamiltonian density have dimensions 



H 
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, n = 


' 1 ■ 




4 


x\ 




X]_ 



Thus all allowed terms without fermion fields are ||GP 92|| 



2 i2 

m , 



Including fermion fields one obtains ||GP 92|| 



_L^{+)tr,^(+)^(+)tr2V;(+) 



(5.19) 
(5.20) 

(5.21) 

(5.22) 

(5.23) 

(5.24) 



(r, Fi and r2 are some Dirac matrices). Unfortunately, this is not the 
whole story. Already free LF-field theory is nonlocal in the longitudinal 
direction [e.g. for scalar fields because of the fundamental commutator 
[0(a;~, x+), 0(?/~, x+)] = —^e{x~ — y~) and for fermions because an inverse 
derivative of 9_ appears in the kinetic energy term (|2.28|) ]. Thus longitudinal 
locality is no longer a restriction on the functional form of the possible terms 
in the Hamiltonian. As a result, any of the operators in the above catalog 
may be multiplied by arbitrary functions of ratios of longitudinal momenta! 
In fact, there are examples known where complicated functions of ratios of 
incoming and outgoing momenta, multiplying a four fermion counterterm, 
are necessary to cancel UV divergences [pH- 
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As a result of this infinitely complicated counterterm structure, it seems 
one looses predictive power and one thus might be forced to abandon LFTD 
as a. fundamental theory of hadrons. It should be emphasized that, even if one 
does have to abandon LFTD as a fundamental theory, it might still have many 
virtues in parton phenomenology. However, there has been recent progress 
towards understanding how to apply renormalization group techniques to 
LFTD that may help restore its predictive power ||WW+ 9^ (for an excellent 
pedagogical review, see Ref. | [Fe 94b|| ). 

Unless one works with the full Hamiltonian, nonperturbative bound state 
calculations in QCD almost inevitably violate gauge invariance |. Therefore, 

' m are excluded. 



®See the discussion in Ref. | GP "9^ why terms with negative powers of ', 
^Lattice gauge theory being the only exception. 
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if one wants to derive a constituent picture from QCD, one is forced to allow 
that explicit gauge invariance is violated: gauge invariance becomes a hidden 
symmetry ||WW+ 9^ , [Pe 94 b|] . In LFTD one introduces cutoffs that violate 
symmetries which normally prevent a constituent picture from arising (gauge 
invariance and full Lorentz invariance). In a sense, the counterterm functions 
that complicate renormalization offer a possible resolution of apparent con- 
tradictions between the constituent picture and QCD | Pe 94b|] . 

The technique to remove cutoff dependence from physical results is renor- 
malization: for example, the functions of momentum fractions that appear 
in the relevant and marginal operators can be fixed by demanding covari- 
ance and gauge invariance in physical observables. An alternative way to 
fix the marginal and relevant counterterms is coupling constant coherence 
(CCC): one insists that functions appearing in non-canonical relevant and 
marginal operators are not independent functions of the cutoff, but depend 
on the cutoff implicitly through their dependence on canonical couplings 
[Fe 94a| , [Fe 94b|| . This automatically fixed they way in which new variable 
evolve with the cutoff, and it also fixes their value at all cutoffs if one insists 
that the new counterterms vanish when the canonical couplings are turned 



of | Pe 94b|| . Remarkably, in the examples studied in Ref. [ PW 93 , Pe 94a 



this procedure provided the precise values for the non-canonical terms that 
were required to restore Lorentz covariance for physical observables. For an 
explicit example for CCC, the reader is referred to Ref. ||Pe 94b| . 



Even without assuming CCC, one can employ renormalization group tech- 
niques [ Wi 75 1 to help determine the counterterm functions: using the pow- 
erful tool of try and error, one makes an ansatz for these functions, which one 
can improve by repeatedly applying renormalization group transformations 
to the effective LP Hamiltonians with these functions included. The fact 
that the QCD Hamiltonian should be an ultraviolet stable fixed point under 
these transformations can be exploited to improve the original ansatz for the 
counterterm functions ||WW+ 9^ - Probably, such transformations alone are 
not sufficient to completely determine the renormalized LP Hamiltonian for 
QCD, but one can improve this approach considerably by using perturbation 
theory, CCC and perhaps phenomenology to guide the ansatz functions used 
in the renormalization group approach to LP Hamiltonians. 

Another promising idea in the context of LPTD is the similarity trans- 
formation ||GW 93|] . Whenever one derives an effective Hamiltonian by elim- 
inating states above a certain cutoff perturbatively, one faces small energy 
denominators, and thus large and uncertain corrections, for states that close 
to (and below) the cutoff. This feature makes it very difficult to repeat- 
edly apply renormalization group transformations because matrix elements 
of states near the cutoff are large. To resolve this problem, Glazek and Wil- 
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son have suggested to apply a cutoff to energy differences instead of to single 
particle energies. By construction, this resolves the problem of small energy 
denominators, but it also provides a band diagonal Hamiltonian. The sim- 
ilarity transformation exploits this type of cutoff and thus provides a way 
to apply renormalization group techniques to LF Hamiltonians (and other 
many body problems) ||GW 93| . 
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Chapter 6 

Summary, Conclusions and 
Outlook 

LF field theory is a very promising approach toward calculating correlation 
functions along a light-like direction. Such correlation functions appear in 
the theoretical analysis of a variety of hard scattering processes, such as 
deep inelastic lepton-hadron scattering and asymptotic form factors. Prob- 
ably the most intriguing and controversial property of LF Hamiltonians is 
the triviality of the ground state. Recent developments indicate that LF 
Hamiltonians must be regarded as effective Hamiltonians in the sense that 
some of the interactions acquire nonperturbative renormalizations with coef- 
ficients proportional to vacuum condensates. So far one understands the LF 
vacuum and is able to construct the effective LF-Hamiltonian only in a few 
toy models. However, in these examples only a finite number of condensates 
are necessary to completely specify the Hamiltonian. It would be extremely 
useful if one could construct and approximately solve such an effective LF 
Hamiltonian for QCD, not only for the analysis of hard processes, but also 
for our understanding of low energy QCD: due to the triviality of the LF 
vacuum, a constituent picture makes sense and an effective LF Hamiltonian 
for QCD would offer the opportunity for deriving a constituent picture as an 
approximation to the QCD bound state problem. 

Three main stream directions can be distinguished in the endeavor to- 
ward constructing a LF Hamiltonian for QCD: First, a fundamental approach 
where all zero-modes are included as dynamical degrees of freedom. Second, 
an effective approach where one attempts to absorb all zero-modes and as- 
sociated vacuum effects into effective interactions and coupling constants. 
Third, the LF Tamm Dancoff approach, where not only vacuum effects but 
also effects from high energy and high Fock components are "integrated out" 
and absorbed into effective interaction terms. In the fundamental approach 
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(this includes all formulations of LF field theory where explicit zero-mode 
degrees of freedom are included) the vacuum as well as the physical particle 
states are complicated and one partly looses the dynamical advantages of the 
LF framework. While it seems easier to construct the Hamiltonian than in 
the other two approaches, the main difficulty of the fundamental approach 
lies in the fact that the equations of motion are extremely complicated. It is 
not clear whether such an approach provides any computational advantage 
over a conventional Hamiltonian approach. Nevertheless, it is very useful to 
pursue this approach further in order to provide a solid theoretical basis for 
other, more practical, approaches to LF field theory. For example, studies 
that include zero-modes can be useful for deriving an ansatz for the effective 
LF Hamiltonian in the large volume limit. 

The LF Tamm DancojJ approach corresponds to the other extreme. The 
vacuum is trivial and the physical particle states are very simple — by con- 
struction they contain only the low energy effective degrees of freedom. A 
major virtue of this approach is that it stays close to physical intuition and 
thus potentially offers a connection between the constituent picture and QCD 
Pe 94b|| . While the LF Tamm Dancoff approach is thus very appealing from 



the intuitive point of view its main disadvantage is the enormous complexity 
of the effective Tamm Dancoff Hamiltonian. In principle an infinite num- 
ber of counterterms are possible. These counterterm functions are heavily 
constrained by imposing Lorentz covariance on physical observables or by 
demanding cancelation of unphysical divergences. However, so far it is not 
clear to what extend one can employ renormalization group techniques to 
constrain the possible interactions to the point where only a few (instead of 
infinitely many) free parameters enter the LF Tamm Dancoff Hamiltonian 
of QCD. The second {effective, in the sense of zero- mode free) approach to- 
ward constructing the LF Hamiltonian for QCD stands in between the other 
two in several respects. The vacuum is trivial but physical particles will in 
general have a complicated wavefunction. Some of the interactions in the 
effective LF Hamiltonian have coefficients proportional to vacuum conden- 
sates. Those can either be regarded as free parameters or (in some cases) 
they can be determined from self-consistency conditions. Surprisingly, in 
those cases where the construction of such an effective Hamiltonian has been 
accomplished, already a finite number of condensates is sufficient to specify 
the Hamiltonian. This is a very encouraging result. Perturbative calcula- 

"'^This approach should not be c onfused with the standard QCD-sum rules approach 
to the strong interactions | NS+ 81 ], where one does not solve a Hamiltonian and where 
practically all the dynamics is buried in the condensates. Hence it is not surprising that 
less condensates are necessary as an input in the LF effective Hamiltonian approach than 
in the sum rule approach. 
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tions up to two loops indicate a similar result for QED, where the two loop 
calculations do not require any countertcrms which are not already present 
at the one loop level. LF perturbation theory in QCD has so far only been 
performed up to one loop. 

Although there has been considerable progress recently, so far none of 
these three approaches has been successful to the point where it was possible 
to construct a useful LF Hamiltonian for QCD. The initial optimism about 
LF quantization, spurred by the very successful application to 1 + 1 dimen- 
sional field theories, was premature. Much work remains to be done before 
LF quantization can be applied to QCD. 

For example, it is still not completely understood to what extend LF 
Hamiltonians, with a trivial vacuum, can account for the phenomenon of 
spontaneous symmetry breaking. The only examples where this subject 
seems to be mostly understood are (j)^ theory in 1 + 1 dimensions and field 
theories in the mean field approximation. It would be interesting to study 
cases where the order parameter for the symmetry breaking does not en- 
ter the Hamiltonian — which is for example the case in the spontaneous 
breakdown of chiral symmetry in QCD. 

A possibly related issue, which requires further study, concerns the non- 
covariant counterterms. In the context of perturbation theory it has been 
shown that a finite number of such counterterms are necessary in the bare 
Hamiltonian to recover full Lorentz covariance for physical observables. How- 
ever, so far it has not been demonstrated that the proposed counterterms are 
sufficient to restore Lorentz covariance for physical observables in a nonper- 
turbative calculation. 

Within the context of LF Tamm-Dancoff it is still necessary to demon- 
strate that the renormalization group, combined with constraints from 
Lorentz invariance, is sufficient to fix the infinite number of counterterms 
which are possible on general grounds. 

For the transverse lattice approach to be useful, it must be shown that 
the fermion sign problem, which usually limits Hamiltonian Monte Carlo cal- 
culations with fermions considerably, is tractable. Since vacuum fluctuations 
are suppressed in LF quantization, any sign problems arising from vacuum 
diagrams are trivially absent. While this is a very encouraging observation, 
it resolves only part of the problem — sign problems arising from exchange 
diagrams within a hadronic state are of course still there. Another difficulty 
for transverse lattice calculations occurs because gauge invariance on such a 
lattice requires the introduction of 1 + 1-dimensional link fields. One must 
learn to work with these "nonlinear sigma model" degrees of freedom in the 
context of LF quantization before one can apply the transverse lattice to 
QCD. 
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Besides QCD oriented applications of the LF formalism, it may turn out 
to be very useful to consider phenomenological and/or more nuclear physics 
oriented applications as well. For example, it may be interesting to reconsider 
the pion contribution to nuclear structure functions |[FP 83|| from the point of 
view of LF quantization. On the one hand, this could be helpful in clarifying 
the role of binding effects in such calculations. On the other hand, such 
works may help to demonstrate the usefulness of LF quantization to people 
who are not directly involved in the field. 

LF quantization is very closely related to the infinite momentum frame 
formulation of field theory. Intuitively one would thus expect that the LF 
formulation of QCD offers a new theoretical approach to relativistic heavy 
ion collisions. So far, this connection has been exploited only very little 
[ |MV 9411 . 
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Appendix A 

The Dirac Bergmann 
Formalism 



In this appendix, a brief introduction into the Dirac- Bergmann quantization 
procedure ||Di 5CI| , Pe 56| , [HR 76|| is given. Quite generally, it replaces the 
canonical quantization procedure in the presence of constraints. However, 
it can also be used to derive the correct fundamental commutation relations 
for theories, where the Lagrangian contains at most linear terms in the time 
derivative. [| 

This is for example the case for many field theories, when expressed in 
terms of LF- variables ||Su 82[| . For example, for a noninteracting massive 
scalar field in 1 + 1 dimensions one obtains 



C = d4 



m 
T 



(A.l) 



which contains no terms quadratic in Naive canonical quantization, i.e 

5C 



with 



yields 



TT 



5d+4 



y 



dx T4. 



_m 
dx — ( 
2 ^ 



-I [P+,9_ 



2m^ 



(A.2) 
(A.3) 

(A.4) 
(A.5) 



^It should be noted that, m the latter case, alternate treatments are possible as well 
IFJ 88|1. 
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Eg. ( |A.5|) differs from the (correct) Euler-Lagrange equation 



ra 



by a factor of two ||CR+ 73| . Tliis mistake arises because 



Clearly, 
-2d+d_ 

the kinetic term in Eg. ( |A.1| ) is only linear in the time derivative 9+. Thus the 
equation relating the canonical momenta to the fields (|A.2|) is a constraint 
equation since it contains no time derivative and therefore the phase space 
variables 7r(x) and (j){x) for a given time are not independent. 

Quantizing a system with constraints is a nontrivial task. Fortunately, the 
Dirac-Bergmann algorithm provides a step by step prescription for the proper 
quantization procedure. The basic steps of this procedure will be illustrated 
in an example below. To keep the discussion simple, zero-modes will be 
deliberately left out in the discussion. A complete discussion, which includes 
zero-modes, can be found in Refs. [|HK+ 89| , |HK+ 91b| , |HK+ 92b| , |HK+ 91a 



Furthermore, the discussion here will be restricted to a system with a finite 
number of degrees of freedom (which can, for example, be obtained from 
Eq. (|A.l|) by discretizing the x~ direction) 



N 



(A.6) 



-Aji (the symmetric part of Aij corresponds to a total time 



with Aij - 

derivative and can be subtracted). The canonical momenta are given by 



0(pi j 



(A.7) 



with Poisson brackets 



{(f)i,TTj} = 5ij 

= {■rc^,^Tj} = 



(A.8) 



Eq. (|A.7|) does not contain any time derivative, i.e. it should be considered 
as a constraint 



Xi = TTi - ^ A 







(A.9) 



i = 1,...,N. Here ~ means weakly vanishing, i.e. a constraint on the 
physical phase space. The canonical Hamiltonian is constructed as usual 



H, = Y^<j>,'K,-L = Y^V{<P, 



(A.IO) 



The constraints ( |A.9| ) have nonvanishing Poisson brackets with H^. 

{X^^H,} = -V\ct>i) (A.ll) 
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as well as among themselves 

{X^, X,} = A, - Aji = 2Aij. (A.12) 

Thus if the time evolution would be generated by Poisson brackets with 
the the theory would be inconsistent because the constraints would not be 
satisfied at all times. To remedy the situation one adds Lagrangian multipli- 
ers to He, yielding the primary Hamiltonian 

Hp = H, + Y.\iXi (A. 13) 

i 

and demands strong vanishing of the Poisson bracket of Hp with the con- 
straints 

= {x^, Hp} = {x^, H,} + J2 ^kBk^ (A.14) 

k 

where Bki = 2Aki. To simplify the discussion, let us suppose that B^^ exists. 
Then one can satisfy Eq. (|A.14|) by choosing 

Xk = J2iB-')k^{x^,H,}. (A.15) 

i 

The primary Hamiltonian thus reads 

Hp = H,-Y. XriB-% {xj, H,} . (A.16) 

Let us now introduce the Dirac brackets between X and Y 

{X, Y}^ = {X, Y}-Y: {X, X.} iB-% {xj, Y} . (A.17) 

By construction one has 

X = {X,H,}^ = {X,Hp}. (A.18) 

Actually, the Dirac bracket of any operator with any of the constraints van- 
ishes identically 

{X,Xk}D = {X,Xk}-J2{X^X.}iB~%{x,,Xk} 

= {X,Xk}-J2{X,X^}{B~%B,k = 0, (A. 19) 

^For the LF-Lagrangian this is actually not the case. There is one zero eigenvalue — 
the infamous zero-mode — which has to be treated separately. The resulting procedure is 
known as the modified Dirac-Bergmann algorithm. 
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i.e. in particular {xi, Xj} = 0- sense, the Dirac brackets take all the phase 
space restrictions from the constraint equations automatically into account. 
This is in sharp contrast to the Poisson brackets, which are calculated as 
if all the TTi's and (pi's were independent {X,Y} = |f |^ - g^fj- It 
thus seems natural to use Dirac brackets, instead of Poisson brackets, when 
identifying classical brackets with quantum commutators 



{X, Y}^ I [X, Y] 



(A.20) 



in the quantization process. A more thorough discussion on this subject 
can for example be found in Ref . |^u 82|| . Here we are more interested in the 
consequences of Eq.( |A.20]) . For this purpose, let us evaluate the fundamental 
Dirac brackets 



]i D 



kl 



kl ^ 



(A.21) 



Roughly speaking, the reduction of the number of independent degrees of 
freedom by a factor of two manifests itself in a factor 1/2 in the Dirac bracket, 
and after applying Eq. (|A.2CI| ) the factor 1/2 also appears in the quantum 
commutator between vTj and ipj. For the LF-quantization of scalar fields this 
implies 



vr X 



(A.22) 



instead of Eq.( |A.3| ). Clearly, this remedies the abovementioned ([A.5|) prob- 
lem with the extra factor of 2 in the LF equation of motion generated by P+ 

From the physics point of view [ |Gr 94a]| , the whole difficulty in quantiza- 
tion with constraints could be avoided if it were possible to choose degrees of 
freedom which are compatible with the constraints. For the above example 
this is actually possible, since the constraint (|A.9|) is linear in the fields. Let 
us thus make the ansatz 



a 



(A.23) 



where a is a constant and the f3ij are determined by requiring a vanishing 
Poisson bracket between 6^ and the constraints 



O = {0f,Xi}=a 



k 



(A.24) 
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i.e. Pij 
i.e. 



{A The normalization a is fixed by demanding that 



a 



a 



+ ^)..J2^jk'^k 



2a(j)i 



(A.25) 

yielding a = \. By construction, (pf does not "see" the constraint and one 
can apply canonical quantization rules directly 



TT,; 



and thus 



(A.26) 



(A.27) 



In the continuum limit this results in Eq. (|A.22|) . 

As already indicated above, there are more points that need to be dis- 
cussed before the Dirac-Bergmann procedure for constructing the LF Hamil- 
tonian is complete. In particular, one has to address the issue of zero-modes. 
In the LF example, the analog of the matrix Aij is the differential opera- 
tor 9_. Fields which are independent of x~ are annihilated by (9_ and thus 
correspond to eigenvectors with eigenvalue zero (on a finite interval, with 
periodic boundary conditions, these zero-modes have to be considered for 
a complete formulation of the theory). In such a situation one first has to 
project on the Hilbert space orthogonal to the zero-modes before the simpli- 
fied procedure above can be applied. The resulting modified Dirac-Bergmann 
procedure is quite involved and has been discussed extensively in the liter- 
ature ||HK+ 89| , |HK+ 91b| , [HK+ 92 b| , |HK+ 91a|| . The basic difficulty arises 
because the constraint equation for the zero-mode is nonlinear. For exam- 
ple, in 0^ theory in a "box" with periodic boundary conditions in the x~- 



direction, integrating the Euler Lagrange equation 
yields f|HK+ 92ai PV 93| , PV 94|| 



-2d^d+(f) = rrrcj) + 



over X 



m 



dx 



+ ^1 dx~ 



0. 



(A.28) 



In the first term in Eq. ([A.28|) only the zero-mode is projected out but in 
the second term higher modes contribute as well. Because the constraint 
equation ( |A.28| ) is nonlinear, the resulting quantum theory is as complicated 
as the formulation in usual coordinates. So far, it is not clear whether any 
dynamical simplifications (like "freezing out" of the zero-mode) arise in the 
infinite volume limit. 
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